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ABSTRACT 


This  is  an  investigation  of  the  field  due  to  a  general  point  source  of 
energy  in  an  isotropic,  elastic  solid  with  a  free  surface.  The  paper  is  divided 
into  three  parts. 

In  Part  I  we  are  concerned  with  the  development  of  new  plane  wave 
representations  for  the  fundamental  solutions  of  elastodynamics.  There  are 
two  types  of  situation  involved;  we  have  the  simpler  type  involved  in  the  case 
of  a  steady  point  source  which  moves  steadily  with  any  constant  velocity  in  an 
elastic  medium,  this  type  involves  superposition  of  plane  waves  with  respect 
to  a  single  parameter,  and  we  have  the  more  complicated  transient  problem  in 
which  a  point  source  is  set  up  at  a  given  moment,  and  thereafter  moves  at 
constant  velocity,  without  change  of  strength. 

In  Part  II  we  make  use  of  the  new  representation  for  the  field  of  a  steadily 
moving  source  in  the  calculation  of  fields  and  displacements  in  the  presence 
of  a  free  surface  and  in  Part  III  we  do  the  same  for  the  transient  source.  We 
discuss  in  some  detail  the  application  of  the  new  approach  to  the  case  of  a 
vertical  load,  to  a  horizontal  load,  and  to  a  couple  of  arbitrary  orientation, 
and  we  give  a  general  discussion  of  the  singularities  to  be  expected  for  the 


general  point  source. 
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THE  USE  OF  SINGULAR  INTEGRALS  IN  WAVE  PROPAGATION  PROBLEMS;  WITH 
APPLICATION  TO  THE  POINT  SOURCE  IN  A  SEMI-INFINITE  ELASTIC  MEDIUM. 


M,  Papadopoulos 


General  introduction; 

In  this  paper  we  shall  consider  the  formation  and  the  propagation  of  tremors 
on  the  plane  surface  of  an  elastic  solid.  Although  considerable  research  has 
been  performed  in  this  subject,  (see  a  complete  bibliography  for  pre-1955  work 
listed  by  Ewing,  Jardetzky  and  Press,  (1957),  and  for  later  work  listed  in  the 
1962  edition  of  Cagnlard),  the  work  to  be  described  here  is  more  general  than 
that  of  Pekeris,  or  of  Cagniard  (1939),  and  it  introduces  a  technique  which  has 
been  found  useful  in  the  study  of  diffraction  (Papadopoulos  1963e). 

With  others,  Pekeris  (1955a,  1955b,  1957,  1958)  has  solved  a  number  of 

problems ,  each  concerned  with  the  setting-up  of  a  point  source  of  energy  by  the 

sudden  application  of  a  constant  force  or  couple.  Both  Pekeris,  and  Cagniard 

(1939),  have  introduced  the  same  kind  of  mathematical  method.  Cagniard,  however 

limited  himself  to  situations  with  an  axis  of  symmetry  normal  to  the  free  surface, 

whereas  Pekeris  (1958)  has  given  one  example,  that  of  the  buried  torque-pulse, 

where  this  property  of  symmetry  is  absent.  In  his  paper  Pekeris  derived  the  form 

for  the  vertical  component  of  the  surface  displacement  without  giving  any  of  the 

detailed  results  for  the  horizontal  components.  Moreover,  although  he  used 
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transform  methods  to  develop  his  solutions,  his  avoidance  of  the  delta  function 
in  the  definition  of  the  point  force  and  in  the  calculation  of  its  associated  field 
introduces  a  note  of  inconsistency  in  the  mind  of  the  reader. 

The  problem  to  be  solved  here  involves,  in  its  most  general  form,  the 
sudden  appearance  of  a  point  source  which  is  moving  with  an  arbitrary  constant 
velocity.  There  is  no  restriction  oii  this  velocity,  beyond  that  implied  in  not 
allowing  this  forced  action  to  cany'  the  source  through  the  surface  without 
specification  of  further  source  behaviour.  The  prototype  of  the  energy  source 
is  a  point  force  of  step  function  time  dependence  and  of  arbitrary  direction . 

Apart  from  implying  that  the  most  general  source  can  be  obtained  from  the  point 
force  by  a  series  of  linear  operations  (vector  addition,  differentiation  and 
integration),  we  have  no  restriction  on  the  type  of  source,  and  we  have  no 
restrictions  on  symmetry  or  orientation. 

If  the  moving  source  is  taken  in  an  infinite  medium  without  a  free  surface, 
the  extention  of  known  fixed  source  solutions  to  describe  the  moving  field  is 
simple.  This  is  because  the  field  is  composed  of  sheer  and  compressional  fields 
which  travel  independently  and  which  may  be  described  in  terms  of  retarded 
potentials  (de  Hoop  1958).  It  is  also  easy  to  take  a  specific  source  and  to  find 
the  moving  field  by  transform  methods  as  described  by  Eason,  Fulton  and  Sneddon 
(1956).  'When  coupling  between  shear  and  compression  effects  is  enforced  at  a 
free  surface,  retarded  potentials  are  certainly  of  no  use,  and  the  most  direct 
approach  seems  to  involve  taking  the  known  field  of  a  transient  source  at  a 
fixed  point,  and  then  superimposing  the  effects  of  a  staggered  distribution  of 
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such  sources  along  the  line  of  motion.  Payton  (1962)  has  used  such  a  method, 
but  when  It  is  applied  to  the  point  source  moving  with  no  acceleration,  it  seems 
to  be  an  overcomplicated  approach.  Physical  Interest  in  the  case  of  a  source 
which  moves  on  a  surface  is  quite  clear;  for  source  motion  inside  a  solid,  there 
is  Interest  not  only  in  possible  earthquake  motions,  but  in  the  effects  of  lines 
of  explosive  charge  due  to  the  finite  propagation  velocity  of  the  explosion,  and 
even  in  the  enveloping  effects  of  a  sequence  of  discrete  explosions  as  used  in 
sequential  rock  blasting. 

The  method  to  be  described  here  will,  in  its  generality,  involve  us  in  the 
calculation  of  known  results  for  the  fixed  source.  In  the  process  we  find  that  the 
results  of  Pekeris  and  Lifson  (19  5T),  for  the  vertical  force  are  incorrect  in  the 
vertical  displacement  component.  The  method  is  of  special  interest  because 
specific  parts  of  the  field  are  picked  out  without  ambiguity;  not  only  can  we 
pick  out  the  singular  parts  of  the  field  for  the  general  point  source,  but  we  can 
pick  out  the  stages  in  the  development  of  the  head  wave  field  as  well. 

There  are  three  stages  in  the  discussion.  In  Part  1,  it  is  shown  that  we  may, 
in  an  Infinite  solid,  represent  the  effect  of  a  point  force  by  an  Integral  superposition 
of  plane  waves.  Both  real  and  complex  plane  waves  have  a  part  in  these 
superpositions.  The  geometrical  envelope  of  the  waves  linked  with  a  specific 
source  forms  a  singular  surface,  or  wave  front,  and  within  this  front  the  field 
at  every  point  is  found  to  be  defined  in  terms  of  values  on  complex  conjugate 
characteristics,  while  the  field  outside  is  similarly  associated  with  real 
characteristics  (or  plane  waves).  In  separate  sections  we  take  the  case  of  the 
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steady  source  which  moves  with  constant  velocity,  this  being  a  steady  state 
problem  in  a  moving  system,  and  the  case  of  the  source  which  appears  and  then 
moves  steadily.  The  first  of  these  leads  to  integral  superpositions  of  plane 
waves  with  respect  to  one  parameter:  it  is  described  in  some  detail  in  order  to 
simplify  the  description  of  the  transient  problem  in  which  two  integration 
parameters  are  needed. 

For  those  who  are  familiar  with  the  problems  of  wave  propagation,  the  value 
of  having  a  plane  wave  representation  for  a  source  field  is  that  the  effect  of  the 
plane  surface  may  be  found  by  treating  each  component  wave  as  being  reflected 
and  refracted  independently  of  the  other  component  wave.  A  formal  superposition 
of  these  reflected  and  refracted  waves  gives  the  total  effect  of  the  free  surface. 

In  Part  2  we  discuss  the  structure  of  the  field  of  the  steadily  moving  source, 
when  it  travels  horizontally,  and  we  calculate  general  details  of  the  surface 
displacement  field.  The  special  case  of  a  fully  supersonic  point  load  moving  on 
the  surface  has  been  described  previously  (Papadopxjulos  1963a). 

In  Part  3  we  examine  the  unsteady  source  field.  We  discuss  in  detail  the 
surface  displacement  for  a  vertical  force,  for  a  horizontal  force,  and  for  a  'double 
force  with  moment*  with  a  specific  orientation,  this  being  regarded  by  seismologists 
(Keilis-Borok  et  al. ,  I960}  as  a  typical  model  for  the  focus  of  a  tectonic 
earthquake.  The  complete  displacement  field  for  an  arbitrary  point  source  may 
not  be  found  without  quadrature,  but  explicit  algebraic  forms  for  all  components 
are  derived  both  for  the  initial  displacements  linked  with  the  arrival  of  primary  P 
and  S  waves,  and  for  those  linked  with  the  airlval  of  head  waves  of  shear,  and 
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for  the  singular  surface  waves  due  Independently  to  the  compressional  and  the 
shear  part  of  the  primary  source  when  on  the  surface.  Of  special  interest  is 
the  discussion  of  the  relative  strength  of  all  the  singularities  which  may  appear. 
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Part  1;  Plane  wave  representations  for  fundamental  solutions  in  elastodynamics. 
Section  1;  The  steadily  moving  source 

The  method  to  be  described  here  was  developed  in  the  study  of  acoustic 
energy  sources  (Papadopoulos  1963b).  A  point  source  of  unit  strength  moves 
along  the  z-axis  with  a  constant  velocity  a;  the  velocity  potential  ^  associated 
with  such  a  singularity  is  the  solution  of  the  equation 


<t>  =  -6(x)  6{y)  6{z  -at)  =  -6(x)  6(y)  6(t-z/a)/a 


{lay 


Since  this  situation  involves  a  steady  motion  in  the  z-direction  we  may  infer  that 
the  velocity  potential  depends  only  on  three  independent  variables,  namely  x,  y 
and  T  =  t  -  z/a  .  Equation  (la)  now  reduces  to  the  form 


<|>  =  -  6(x)  6(y)-6(T)/a  , 


(lb) 


2  2 

with  Y  =  3c(a  -c  )  ^  ;  in  equation  (lb)  we  can  see  the  well-known  change 
from  hyperbolic  to  elliptic  form  as  we  permit  a  to  pass  from  supersonic  values, 
with  a>  c,  to  subsonic  values  with  a  <  c  .  The  equations  (1)  are  appropriate 
for  any  real  value  of  a  .  When  a  s  0  we  have  a  static  situation,  and  when 
a  =  eo,  we  have,  after  multiplying  the  right  hand  side  by  a  ,  the  equation 
defining  the  potential  of  an  infinite  line  source  of  unit  strength  per  unit  length. 

For  supersonic  values  of  a,  the  transverse  velocity  of  propagation,  y> 
real,  and  the  solution  of  equation  (lb)  has  been  given  in  the  form 
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(j,  =  -  Rl  /  - - ^-r -  dp  .  (Za) 

4tt  a  -00  (1-p  )^[yt-px-(1-p  y  sgny] 

For  subsonic  values  of  a,  the  transverse  velocity  is  imaginary  (y  =  ijx,  say,) 
and  the  solution  of  equation  (lb)  has  been  given  in  a  slightly  different  form 
with 

<t)  =  -^Rl/  - ^-T - ^ •  (2b) 

4iT  a  -ioo  (1-p  )2  [i^iT  -  xp-ysgny(l-p  )^] 

The  main  value  of  these  representations  is  that  away  from  the  source,  the  wave 
nature  of  the  associated  field  is  evident.  Corresponding  expressions  associated 
with  strain  nuclei!  will  now  be  considered. 

The  eventual  problem  to  be  examined  involves  the  setting  up  of  elastic 
disturbances  in  a  homogeneous  isotropic  medium,  of  density  p  and  with  Lam^ 
constants  X.  and  p.,  in  a  half  space  0  with  the  plane  y  =  0  a  free  surface. 
We  set  up  the  velocity  vector  v  in  terms  of  velocity  potentials  ^  and 
(=  M  +  Bj_  +  Ck  ),  such  that 

V  =  V(^  +  Vxip  ,  (3a) 

and 

V*  iji  =  0  .  (3b) 

Given  a  body  force  F,  the  stress -strain  relations  and  the  equations  of  motion 


reduce  to  the  forms 
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and 


2  w2 

c,  V 


2  v72 


^2  ,  9  — 

^  i  =  iT  -  ’ 


(4) 


2  2 

with  pCj  =  \  +  2  (i  ,  pc^  =  K  >  and  being  the  velocities  of  propagation 
of  P  and  S  waves  respectively.  When  a  force  p(}Q  +  Xi  +  Zk )  is  applied 
steadily  to  the  moving  point  x  =  y  =  0,  z<_at,  F  has  the  form  given  by  the 
equation 


F  =  (la  +  YJ.  +  Zk )  6  (x)  6  (y)  6  (z-at)  , 

=  ()a  +  Y1  +  Zk )  6  (X)  6  (y)  8  (T)/a  ,  (5) 


and  the  equations  (4)  reduce  to  the  transverse  forms 
2 

°1 


and 


’  1 

8^ 

8^ 

r  2  2  21 

9  .-C.  1  9 

2 

L^i 

8r^ 

'ax^' 

8y^ 

^2'^^2‘^  2^2 
8x  8y  a  8t  _ 

♦=~V-r,  (6a) 


i  .ii  £  z  1  r 

2^2  8x  ^2 

Y2  8y  J  L 

with  =  aCjCa^-Cj)"^,  and  Y2  =  (a^  -  c^ ) 


z  z  z 
-lit  -iitJ— li 

8x" 


i=-^VxE,  (6b) 


Equation  (6a)  has  a  simple  solution,  given  as  the  Inverse  of  a  triple  Laplace 
transform  by  the  equation 
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but  because  the  integrand  is  homogeneous  in  the  three  transform  variab'.es  we 
may  put  V  =  -  sp/y^,  p  =  -sq/y,  so  as  to  reduce  equation  (7)  to  the  form 


Yj  (pX+qY+LZ)  exp  s[Y,T-px-qyJ/Y, 

<!>  =  ■  I  I  J  Jj  - r - - - - - r - - - i  dpdq  , 

aCj  (Ziri)  -ioo  -ioo  (1-p  -q)(p  +q  +L) 


with  L  =  Y|/^  •  The  integration  with  respect  to  s  may  now  be  carried  out, 
so  that 


<}>  = 


c^  a{2iri) 


(pX-!-qY+  LZ)  6[YjT-px-qy]/Y^ 

JJ  _  _  __  2  2 

-00  (1-p  -q  )(p  +q  +L  ) 


dpdq  ; 


(8) 


the  presence  of  the  delta  function  is  meaningful  only  as  long  as  the  function 
(YjT  -  px-  qy)/Yj  is  real.  In  the  work  being  considered  it  is  not  always  possible 
to  assume  this,  and  the  delta  function  must  be  regarded  as  having  an  alternative 
singular  form,  with  (8)  taking  the  form 


^  = 


'1 


00  00 


^  2  2 
4ir  c^a 


«!/  / 


1  (PX  +  qY  +  LZl 


2  2  2  2  2 

«  irCl-p  -q  )(p  +q  +L  )(Y2T-px-qy) 


dpdq 


(9) 


Away  from  the  point  singularity,  <})  satisfies  the  wave  equation,  and  the 

2  2 

integration  parameters  must  be  restricted  to  the  circle  p  +  q  =  1  .  Thus  for 
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,  ^  ,  ,,,  2  4  and  in  turn  this 

arbitrary  p,  q  may  take  only  the  values  q=*(l“P) 

means  that  the  q-lntegration  involves  only  the  calculation  of  residues  at  these 

2  ^ 

points.  We  specify  the  branch  of  (1-p  by  taking  the  p-integratlon  path 
along  the  real  p-axls,  with  indentations  taking  it  above  the  point  p  =  1  and 
below  the  point  p  =  -1  .  With  this  convention  we  have 


A  =  ^  Ri  (°  [pX-t-d-p^}^  sgn  y  Y  -I- 

“  ”  2  J  •>  i 


LZ 


^  2 
4tt  a 


-00 


2  2 
(l-p  [YjT-px-ysgny(l-p  )^] 


dp 


(10a) 


or,  following  a  change  in  scale  of  the  integration  variable 


-L-  Rj  f  ^  ■  pXt(MSp^l^.Y?gny,t  LMZ 

4ir^a  -00  (M^-p^)^  [Y2T  -  px  -(M^-p^)^  y  sgny] 


dp 


(lla) 


with  M  =  Y2/Yj'^  1  •  These  representations  are  only  correct  for  fully  supersonic 
source  motion,  with  a  >  c^  .  For  the  range  a  <  c^,  the  transverse  velocity 
Y^  is  imaginary,  with  say,  and  the  results  of  simplifying  the  Integral 

(7)  are  slightly  different,  taking  the  integration  paths  for  both  p  and  q  along 
the  imaginary  axis  instead  of  the  real  axis.  To  correspond  to  the  formula  (10a) 
we  have  the  Integral 

♦  =  ItPX t(i-p¥ V;giiyt tfl  . 

4ir  a  -ioo  (l-p  )2  (l,i^T-px-(l-p  )*  ysgnyj 


(10b) 
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The  integral  (11a)  keeps  Its  form  when  a>c^,  but  it  changes  to 

^  Rl/*  trpXHM^-p^Ysqny^t^Uiel - 

4ir^a  -i«0  (M*^-?  (i|i2‘''"P^“(^  -pji^ysgny] 

in  the  case  of  fully  subsonic  source  motion,  with  a  <  c^;  when  y^  =  ^^2  * 

The  three  cartesian  components  of  »J)  satisfy  equation  (6b)  and  may  be 
written  in  the  form 

1  r“  fLMY-Z(l-p^)^  sgn  y]l.  +  [pZ-LMX]L+ [X(l-p^)^ sgny - pY]k 

±  R1 J  - ^ ' - 2  i  — - 

4Tr  a  -00  (1-p  )^[Y2T-PX- (1-P-)*  y  sgn  y] 

(12a) 

when  a  >  0^,  or 

,  _1_  p,  i{[LMY-Z(l  -  p^)^  sgny]i.+  [pZ - LMX]l,-h  [X(l  -  p^)^  sgny -  pY]k} 

Y*  2  ^  y 

4iT  a  -ioo  (1-p  )*  (ip2‘'’“P^”0"P  )*ysgny] 

(12b) 

when  a  <  , 

Notice  here  that  as  the  source  velocity  is  varied  and  the  transverse  velocities 
take  positive  imaginary  values,  the  other  velocity  parameters,  L  and  M,  may 
also  take  imaginary  values. 

The  potentials  which  correspond  to  other  steadily  moving  strain  nuclei!  may 
easily  be  constructed  by  superimposing  the  results  for  simple  forces.  Thus  two 
equal  forces  acting  to  oppose  each  other  in  the  x-dlrection  but  with  a  moment  about 
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the  z-axis  may  be  taken  in  a  limiting  form  as  a  'double  force  with  moment';  the 
potentials  in  this  case  are 

♦  =  - Vn - ‘‘p 

4n  -00  [y2T-PX-(M  -p  )*  ysgny] 

and 

i  =  «  “  /”  dp  . 

4ir  av^  -®®  [Y2'r-Px-(1-P  )^ysgny] 


Similarly  a  combination  of  three  equal  orthogonal  double  forces  without  moment 
gives  the  potentials 


«!>  = 


1 


^  2 


00 


-00  {M^-p^)nv. 


M^.C1  +  L^). 


2  2  - 

T-px-(M  -p)^ysgny3 


dp 


and  ^  =  0  ,  This  being  a  constant  real  multiple  of  the  time  derivative  of  the 
potential  for  an  acoustic  source,  we  might  concentrate  attention  on  this  acoustic 
source  as  being  a  simple  case  which  does  not  produce  shear,  while  to  take 
account  of  shear  the  simplest  case  is  that  of  the  simple  force.  A  detailed 
examination  of  other  strain  nucleii  (e.g. ,  as  listed  by  Love  (1927)  in  the  study 
of  elastostatics ,  or  by  Keilis-Borok  (I960)  in  the  study  of  seismology)  is  not 
needed  here.  Thus  in  addition  to  the  potentials  (11)  and  (12)  we  have  to  refer  to 
the  potentials  of  a  moving  source  of  unit  strength,  in  the  form 

^  =  -- ~Ri  /  ■  -2-  1 1 - ^  r-ji - 

4ir  a  -<»  (M  -p  )*  [y2T  -  p\ -(M  -p)*  ysgny] 


dp  (13a) 
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when  a>  c^,  or 

^  ^  ?1  J  22“  Z  Z  ^  dp.  (13b) 

4ir  a  -00  (M  -p  )*  [ijA^T  “px-(M  p  ysgny] 

« 

Each  of  the  integrals  given  may  be  reduced  to  an  explicit  formula.  The 
process  involves  shifting  the  formal  integration  path  to  the  curve  on  which  the 
wave  function  of  the  denominator  is  real.  On  this  path  the  integral  is  singular, 
the  Cauchy  principal  value  is  either  zero  or  Imaginary,  and  the  only  real 
contribution  to  the  potentials  comes  from  the  residue  at  points  for  which  the 
wave  function  vanishes.  For  supersonic  motion  there  are  two  conjugate  zeros, 
while  for  the  subsonic  case  there  is  only  one.  Real  zeros  of  the  wave  function  do 
not,  in  the  integrals  given,  provide  real  residue  contributions  to  the  potentials. 
Each  one  of  the  integrals  has  a  singular  geometry.  In  all  cases,  the  source 
point  is  a  singular  point  of  the  integral,  and  in  the  supersonic  cases  there  is  a 
singular  (conical)  surface  to  contain  the  source  field.  These  singularities  are 
determined  by  the  simultaneous  vanishing  of  the  denominators  of  the  integrals 
and  their  derivatives  with  respect  to  p  .  The  conical  surface  is  seen  to  be 
the  envelope  of  plane  waves  generated  by  the  source.  The  property  that  complex 
zeros  of  the  wave  function  are  the  only  ones  to  provide  residue  terms  that  matter 
for  the  source  fields  is  equivalent  to  saying  that  these  fields  are  determined  by 
complex  plane  characteristic  surfaces  passing  through  each  point  of  observation. 
The  details  are  similar  to  those  given  previously  (Papadopoulos  1963b). 
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Eacii  plane  wave  in  one  of  the  given  singular  integrals  is  reflected  and 
refracted  at  a  plane  boundary  just  as  if  it  were  a  real  plane  wave.  In  order  to 
fix  ideas  we  shall  consider  the  case  of  a  point  source  moving  at  a  uniform  depth 
h  below  the  free  surface  of  an  elastic  solid  in  the  plane  y  =  0  .  Away  from 
the  source  point  the  potentials  must  satisfy  the  equations 


and 


-  — 


cf  at^ 


•t*  =  0  , 


4)  =  0  , 


(14a) 


{14b) 


for  y  ^  0  ,  while  the  conditions  of  vanishing  stress  at  the  surface  are  that 

at”  ~  4  [<j>  +A  — C  (15a) 

at  yy  T  r  iTyy  2y  xy*  ’  '  ' 

~T  =u[2i  +A  -C  +C  -B  1  =  0,  (15b) 

at  yx  ^yx  xz  XX  yy  zy*  ’ 

~T  =iJi[2«fc  +A  -A  -C  +B  1  =  0, 

at  yz  ^yz  zz  yy  zx  xy*  ’ 


for  y  =  0  .  The  third  condition  may  be  combined  with  (15b}, (14b)  and  the 
wave  equation  to  give  the  simpler  form 


=  0 


(15c) 
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which  implies  that  the  potential  B  Is  an  even  function  of  y  about  the  plane 


y  =  0  . 


In  the  absence  of  boundaries  the  source  field  Is  of  two  slightly  different 
forms.  For  y  >  -  h,  we  have  a  general  source  field  in  the  form 

^  P  ^  A1  +  B  j_  +  Ck 

<j)=  Rl  J — ; — rr - y — 71 - dp,  ±  =R1  j - 71  71 

(M‘^-p‘^)^[Y2T-PX-{M'^-P  )*(y+hj]  (1-p  )*[y2‘^"PX-{1~P  )^(y+l')] 


dp 


(16a) 


with 


pA  +  LMC  +(l-pfB  =  0 


For  y  <  -h,  there  is  the  corresponding  form 

3|c  a|e 

.  F*  /-Ai+BL+Ck 

^  =R1  J—T— [ - 7-YI -  dp,  vjixRl  J - r-T - 71 - dp 

(M  -p  )*[y2T-PX  +  (M  -p  )Hy+hl]  (1-p  }*[y2‘^"P*H1“P  )Hy+h)I 

(16b) 


with 


Id  *  2. 4  Id 

pA  +  LMC  -(1-p  B 


«  0  . 


«  *  Id  * 

It  is  presumed  that  the  functions  F,  A,  B,  C,  F  ,  A  ,  B  and  C  are  real 
functions  of  p  when  the  transverse  velocity  parameters  are  real.  The  integrals 
(11),  (12)  and  (13)  are  those  we  have  specially  in  mind  to  provide  explicit  forms 
for  the  nunerators  of  the  integrals  (16),  but  there  is  no  restriction  to  these.  We 
may  also  consider  derivatives  of  these  integrals  with  respect  to  t,  just  as  we 
may  consider  more  general  linear  operations  with  respect  to  t  as  may  be  needed 
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to  provide  results  for  sources  which  vary  with  time.  From  the  conditions  at 
y  =  0,  it  follows  that  we  may  write  down  the  scattered  field  in  y  <  0  in  the 


form 
<t.  =  Rl/- 


(M^-p^)^p) 


.2  2.4, 


2  2.4 


2  - 


Y,T-px+(M  -p  )2(y-h}  Y  T-px  +  (M  -p  )*y-(l-p  )*h 


f  dp  , 


^=Rl/ 


1 


(iV)^R(p) 


ApL+BpL+Cjjc 


y  T-px+Cl-p*')^  y  -hCM^-p*")^  Y7T-px  +  (l-p^)2  (y  - h) 


2  i 


.2  2.4 


dp 


(17) 

where 

R(p)  =  (1  -  -  2p^)^  +  4(1  -  p^)^(p^  +  L^M^)  ,  (18a) 


the  scattering  coefficients  linked  with  the  primary  P  field  are  defined  by  the 
equations 

Fp  +  FR(p)=  8C1-pVcP^  +  L^M^)  F  , 

Cp  =  -  pAp/LM  =  -4p(l  -  p^)*  (1  -  -  2p^) 

and 

Bp  =  0  , 

while  the  scattering  coefficients  linked  Independently  with  the  primary  S  field 
are  defined  by  the  equations 


I 


(18b)- 
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■17“ 


F-  =  4(M^-p^}Hl-L^M^-2p^)  (pC-LMA)  , 
o 

Bg=BR(p)=  -R(P)  (LMC  +  pA)(l“P^)“*  ,  ► 

[Cg  +  CR(p)]  =  -p[Ag  +  AR(p}]/LM  =  8p(M^-P^)^(l-p^)UpC-LMA)  . 

{18c) 

The  integrals  (17)  represent  the  extra  field  which  has  to  be  added  to  the  primary 
field  (16a)  in  the  region  0>  y  >  -  h,  and  to  the  primary  field  (16b)  in  the  region 
y  <  -h  .  The  field  expressions  for  a  sotarce  of  the  same  form  travelling  on  the 
free  surface  are  to  be  found  by  taking  the  results  in  the  region  y  <  -  h  in  the 
limit  as  h  -*  0  . 


A  general  discussion  of  these  results  will  be  given  in  Part  2. 
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Section  2;  The  transient  point  source. 

We  have  dealt  in  Section  1  with  the  fields  of  steadily  moving  sources; 
these  fields  have  been  derived  in  the  form  of  integrals  of  singular  plane  waves 
with  respect  to  a  single  parameter.  To  give  similar  forms  for  transient  point 
source  fields  involves  us  in  integrals  of  singular  plane  waves  with  respect  to 
two  parameters. 

We  shall  consider  first  the  solution  of  the  equations  (4)  when  the  body 
force  F  is  given  by 

F  =  (XL+ Y1+ Zk)6(x)  6(y)  6(z)  U(t)  .  (19) 

The  equation  for  the  scalar  potential  has  the  solution 


<!>  = 


ioo 

M- 


exoFstt  X-X  -huy+vzKX.X-t-uY+vZ) 


4  jjjj  2  2,,  2,  2  2,,r.2  2  2, 

(2iri)  -ioo  [s  -c^(\  +(i  +v  )][\  +  +v] 


ds  dX  d|j.dw 


(20) 


this  being  the  inversion  of  a  quadruple  Laplace  transform.  The  integrand  is  a 
homogeneous  function  of  the  four  transform  variables,  and  we  may  therefore  put 
X  =  -s<^/Cj,  |x  =  -s  01  y  =  -  s  to  derive  the  expression 


4»  = 


(27ri) 


A  Iff 


(o^X+a^Y  +  OfjZ) 


ioo 


-00 


da^dof^da^  /  expsfc^t-c^x-ft^y-Q^zJ/Cj 


-loo 


00 


■mir 


(«lX+a2Y+Q'3Z)  Cj 


3  2  2  2  2  2  2  2 
(Zir)  -00  Tr(l-«j-ft2-«3)(o'j+cif2+«3)(Cjt-<i'^x-«2y-«3z) 


d®!  dar^  da3 
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This  expression  satisfies  the  wave  equation  away  from  the  origin,  and  for 

2  2  - 

given  a,,  a  is  restricted  to  one  of  the  values  *p=±(l-a;.  -a  .  Thus 
the  whole  integration  with  respect  to  niust  involve  only  residues  at  these 
two  points.  Here  we  define  the  branch  of  the  radical  so  that  the  real  a^-axls 
passes,  with  indentations  as  required,  above  the  point  a^  =  P,  and  below  the 
point  a^  =  “P  .  These  poles  give  the  residue  contribution 


+  a^+  a^Z)  sgn  y 
— f — ; - r  da,  da, 

“ztW’VVJ  *  ^ 


(21a) 


with  «2  “  P  dividing  both  integration  variables  by  m  =  c^/c^. 


^  00  i(ajX+ a^Yf  a^Z)  sgn  y 


(21b) 


2  2  2  - 

with  a2=asgny,  a  =  [m 


In  this  integral  it  is  required  that  the  'slowness'  variables  be  represented  by 
a  point  on  a  spherical  slowness  surface  of  radius  m  .  The  integral  is  therefore 
Invariant  under  rotation  of  these  variables,  this  rotation  may  be  chosen  to 
simplify  further  calculations. 

The  equation  (4) ,  which  defines  the  vector  potential  ij; ,  may  be  solved  in 
a  Eimllcur  fashion  to  give 


Stt  c^  -00 


isgnyKYa^-Zp^)!  +(Za^-Xa^)l+(Xp^-Ya^)k] 


dai  da^ 


(22) 
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with  (3,  =  P  sgn  y  .  Note  that  the  wave  functions  of  both  the  Integrals  (21b) 
and  (22)  are  identical  in  their  behaviour  on  tlie  plane  y  =  0  .  This  property 
has  been  chosen  deliberately  to  simplify  the  matching  of  P  and  S  waves  at 
the  free  surface. 

An  expression  of  similar  form  for  the  scalar  potential  of  a  dilatation  source 
of  unit  strength  is 


<b  =  — R1  fj  — r - - :  da  da  . 

Stt^c  aic^t-a^x-ti^z-ay  sgn  y)  1  3 


Simple  extensions  of  these  formulae  may  be  obtained  to  cover  the  case  of  a 
transient  source  which  moves  with  constant  velocity.  In  this  situation  the  body 
force  is 

F  =  (Xi_+ Yj_+ Z]c)6(x-Vjt)  6(y-V2t)6(z- v^t)  U(t)  .  (24) 


The  effect  of  the  motion  is  to  introduce  into  the  integrand  (20)  an  extra 
homogeneous  factor  s(s+ Xv,  +  (iv_  +  vv  )  ^  .  The  details  of  the  subsequent 
evaluation  are  not  changed;  the  motion  of  the  source  leads  us  to  include  in 
the  integrals  (21b)  and  (23)  and  extra  factor 

Vi  =  OjCOj  -  -  .jVj  -  .jVj)"'  {25a) 


(25b) 
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With  Integration  paths  already  chosen  as  curves  on  which  the  wave  function  is 
real  the  general  set  of  potentials  associated  with  the  appearance  of  a  source 
at  the  initial  point  x  =  z  =  0,  y  =  -h  may  be  written  in  the  form 

iF  Vj 

^[c^t  -  ofjX  -  a^z  -  a(y  +  h)]  *^*1  ^*3  ’ 
i(Aj^+  Bi_+Ck]  V2 

p[c2t  -  flfjX  -  «^z  -  P(y + h}]  *^“1  (26a) 


for  y  +  h>  0,  and 


**  V 1 

'P  ~  ^  If  c^c^t  -  flTjX  -  «^z  +  a(y + h)]  ^*1  ^*3 

l[A\+B*l+C*kJV* 

If  PEc^t  -  «jX  -  OjX  -  a^z  +  p{y  +  h)]  ^*1  ^*3 


(26b} 


for  y  -I-  h  <  0,  the  two  representations  being  needed  to  eliminate  in  the 
numerators  the  function  sgn(y-l-h)  .  (These  expressions,  of  course,  represent 
only  potentials  of  degree  -1,  but  to  give  a  general  form  for  a  potential  of  degree 
-n  we  need  only  differentiate  the  Integrals  (26)  (n  -l)-times  with  respect  to 
t  ). 

Given  a  primary  source  with  potentials  of  this  form,  the  effect  of  a  free 
surface  in  the  plane  y  *  0  is  to  introduce  scattered  potentials,  in  y  <  0  is  to 
introduce  scattered  potentials,  in  y  <  0,  given  by 
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<i*=Rij(r 


c^t  -  OfjX  -  ofjZ  +  o(y  -  h)  *  c^t  -  a^x  -  +  cy  -  Ph  j  ^*3 


and 


‘k=«iir^l 


A  i  +  B  J_  +  C  k 

PT  V~  ..iT. 


A  l+B  1  +  C  k 
s—  S'*"  s  ~ 


®2^  “  “  “32  +  Py  -  ah  ^  c^t  -  a^x  -  +  P(y  "  h)  j  *^“1  *^*3* 


Then,  with  the  function 


„  r.  ,  2  ,  2,2  ^  2  2- 

R  =  (1  -  2aj  -  2«jJ  +  4  ap{aj  +  a^)  » 


(28a) 


appearing  as  a  determinant  of  the  boundary  conditions  (15),  we  find  the  scattered 
amplitudes  of  each  individual  plane  wave  from  the  equations 

Fp  +  FRVj  =:8aP(«j+a3)FVj 
Cp  =  -Apaj/aj  =  -4  Paj(l -2«J -2«^)  F  Vj 

Bp  =  0  (28b) 

and 

Cg  +  CRV^  =  -  (Ag  +  ARV2J/«3  -  A^^J  V2  , 

Fg  s  4(1  -  2afj  -  2ag  )  a(Caj  -  A«g] V2  , 

Bg=BRV2»-pR(A«j  +  C«3JV2  . 


(28c) 
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The  expressions  (28b}  and  (28c)  are  quite  independent,  the  former  represent 
reflection  and  refraction  coefficients  linked  with  an  incident  P  wave,  and  the 
latter  are  linked  similarly  with  an  incident  S  wave* 

In  the  region  0  >  y  >  -h,  the  total  field  is  given  as  the  sum  of  the 
expressions  (26a}  and  (27);  in  the  region  y  <  -h  it  is  likewise  given  as  the 
sum  of  the  expressions  (26b}  and  (27).  The  case  when  the  source  is  formed 
at  the  surface  is  obtained  by  taking  the  field  in  y  <  ~h  and  then  allowing  h 
to  vanish. 

One  restriction  remains  to  be  mentioned  here.  Positive  values  of  v^  will 
bring  the  source  to  the  surface  when  0.  <  t  <  h/v2  .  The  associated  field 
expressions  are  therefore  valid  only  in  a  restricted  region,  and  a  complete 
description  of  the  subsequent  field  has  to  follow  a  choice  of  the  subsequent 


source  motion. 
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Part  2;  The  steadily  moving  source  in  an  elastic  solid  with  a  free  plane  surface. 
Section  1;  The  simplification  of  integral  solutions,  and  the  field  structure. 

We  shall  first  examine  the  means  for  simplifying  integrals  given  in  general 
form  as  in  equations  (1.16)  and  (1.17),  for  the  case  of  fully  supersonic  source 
motion.  For  the  primary  source  potentials  of  equation  (1. 16),  the  evaluation 
Is  performed  by  shifting  the  integration  path  from  the  real  p-axis  to  the  curve 
on  which  the  wave  function  of  the  denominator  is  real.  Thus,  for  the  primary 
dilatation  field  of  equation  1. 16a,  the  wave  function 

Y2T  -  px  -(y  +  h)(M^-p^)*  (1) 

is  real  on  the  hyperbola  p  =  M  cosh(w+i0)  for  -«<  w<  «,  with  6sarctan(y-l-h)/x  . 
2  Z  ^ 

The  branch  of  (M  -p  )*  chosen  has  a  positive  real  part,  with 
2  2  - 

(M  -  p  )^  =  -IM  sinh(w+ 10),  in  order  that  the  wave  function  (1)  be  part  of  a 
distudsance  which  moves  in  the  direction  of  Increasing  y  . 

With  0  chosen  and  the  branch  of  the  radical  defined  in  this  manner,  it  is 
clear  that  the  wave  function  (1)  does  not  vanish  when  r  <  0  .  With  given 
positive  T,  it  has  complex  conjugate  zeros  on  the  hyperbola  paMcosh(MN-iO] 
for  0  <  [x^  +  (y+h)^]^  <  YjT,  and  it  has  real  zeros  if  (x^+(y+h)^J^>YjT  .  It 
follows,  having  shifted  the  integration  path  from  the  real  p-axis  to  the  locus  of 
complex  zeros,  that  there  are  three  distinct  contributions  to  consider.  First  the 
shift  in  path  Involves  us  with  residue  contributions  from  real  poles;  for  the 
primary  source  terms  these  residues  are  strictly  imaginary,  and  do  not  contribute 


to  the  field.  Second,  the  line  integral  along  the  hyperbola  Is  singular  (the  path 
is  shown  In  figure  la),  but  the  Cauchy  principal  value  has  a  zero  real  part.  This 
is  because  the  integration  path  is  symmetrical  about  the  real  p-axls,  and  the 
contribution  to  the  real  part  from  the  lower  half  of  the  path  annuls  that  from  the 
upper  part.  The  third  contribution  is  the  residue  contribution  from  the  complex 
poles;  for  the  primary  dilatation  field  these  are  at  the  points 

y/x  *  i(y  +  h){y2T^  -  M^(x^  +  (y  +  h)^]}* 

P  3  ^ 

X  +  (y  +  h)^ 

2  2  - 

and  these  appear  only  when  YiT  >  [x  +(y  +  h)  J*  >  0  .  The  bounding  surface 

2  - 

for  this  residue  contribution,  a  section  of  the  cone  =  [x  +  (y  +  h)^  >0,  is  a 
singular  surface  of  the  dilatation  field.  It  is  associated  with  double  zeros  of  the 
wave  function,  that  is,  with  points  for  which  both  the  wave  function  and  its 
derivative  with  respect  to  p  vanish;  this  is  the  precise  definition  of  the 
envelope  of  the  individual  plane  waves.  It  marks  the  boundary  surface  for  the 
contribution  of  residues  at  complex  poles,  1.  e. ,  for  the  effects  of  the  conjugate 
complex  characteristics  of  the  wave  equation. 

The  primary  potentials  (1. 16a)  and  (1. 16b)  reduce  in  this  manner  to  the  form 


yy  -  MV  +  (y  +  h)^] 


<|>/2tr  s 


f 

[J^(Pj)]  for  y  +  h  >  0 
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with  p^[x^  +  (y+h)^]^  =  Yj;TX  -  i(y  +  h)  {Y2T^-M^[x^  +  (y  +  h)^]}^sgn(y  +  h)  , 

and  YjT  =  Y2'''/^  ^  ^  ®  > 


while 


2  2  2  ,  _.2 
y,t  -(y+h) 


2  ±jZ-n  =  J 


R1[^p1L+ B(P)L+ C(p)iL]  _  fory  +  h>0 

P-P2 

R1  [A*(p)  L+  b*(pU+  c*{p)  k  ] 


p=p, 


for  y  +  h  <  0  , 
(2b) 


with 


P2[x^  +  (y  +  i^)^]^  =Y2’’x-i(y  +  h)[Y2T^-x^-(y+h)^]2  sgn(y+  h) 


and 


Y2T  >  [x^  +  (y  +  h)^]2  >  0  . 


For  completeness,  it  should  also  be  mentioned  here  that  the  shift  of 
integration  path  should  involve  contributions  from  the  circle  at  Infinity.  These 
make  no  contribution  to  the  field  unless  specific  points  in  the  vicinity  of  the 
source  are  being  examined  or  if  t  -►  0  , 

The  integrals  (1. 17)  are  a  little  more  difficult  to  evaluate.  They  are 
distinguished  by  the  presence  of  the  function 

R(P)  =(1-lV-2pV  +  4(1-p^)^(M^-pV(P^+ •  (3) 

This  is  a  characteristic  function  defining  the  transverse  velocity  y  of  Rayleigh 

lx 

2  2-4 

waves.  It  has  zeros  at  the  points  p=*R  =  y,/Yd>  where  Yu  =  ®C-(a  - c^)  * 

^  R  R  R  R* 
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this  result  being  easily  deduced  from  the  limiting  two  dimensional  situation 
when,  with  a  -♦  <»  and  m  =  c^/c^  , 

R(p)  -  (l-2p^)^  +  4p^(l-p^)^(m^-p^)2  =  0 

is  the  equation  with  roots  p  =  *  actual  velocity  Vj^  of 

Rayleigh  waves. 

The  zeros  of  R(p)  are  poles  of  the  integrals  (1. 17)  ,  They  may  only  be 
linked  with  singular  surfaces  of  the  scattered  field  if  they  are  actually  poles  of 
second,  or  higher,  order,  and  this  is  the  case  only  if  the  wave  function  of  the 
denominator  vanishes  with  R(p)  .  It  is  clear  from  the  form  of  equation  (3) 
that  R(p)  will  not  vanish  except  for  real  values  of  p  with  ipl>  1  .  Other 
(complex)  zeros  which  are  found  after  R(p)  is  made  rational  are  not  zeros  which 
can  be  associated  with  the  branches  defined  for  the  radicals.  The  wave  functions 
for  (1.17)  may  not  vanish,  however,  except  in  the  single  case  with  h  =  y  =  0  , 
y  being  restricted  to  non-positive  values  in  these  expressions.  Singular 
Rayleigh  waves  do  not  occur  except,  on  the  free  surface,  for  sources  set  up  on 
the  surface. 

Notice  that  the  integrals  (1. 17)  are  of  two  types.  The  first  integral  for  ^ 
and  the  second  Integral  for  ({i  are  simple  in  structure,  being  integrals  of  plane 
waves  originating  at  the  singular  image  point  y  =  h,  x  =  0,  z  =  at  .  These 
integrals  are,  in  turn,  reflected  P  and  S  fields.  They  are  evaluated  in  the 
same  manner  as  the  primary  field  expressions,  by  shifting  the  path  of  integration 
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to  the  curve  on  which  the  wave  function  is  real.  For  the  P  field  there  is  nothing 

new  to  discuss;  we  have  only  residue  contributions  from  conjugate  complex 

poles  to  evaluate  and  this  defines  a  field  contained  between  the  surface  y  =  0 

2  2 

and  the  image  cone  Yj"’’  +(y“h)  ]*  for  t  >  0  .  The  reflected  S  field 

2  2  - 

contained  within  the  cone  =  [x  +  (y  -  h)  ]*  is  obtained  from  the  complex 
zeros  of  the  wave  function,  but  there  is  an  extra  complication  to  be  examined. 
The  function  R(p)  contains  branch  points  at  p  =  *M,  with  M  <  1;  the  locus 
of  complex  zeros  of  the  wave  function 

Y2‘*'"P^  +  (4) 

is  bound  to  cross  the  real  p-axis  at  points  with  ipl  <  1  and  [pi  >  M  if 

arccos  M  >  larctan(y  -  h)/x|  >  0  . 

Then,  as  shown  in  figure  (lb),  the  shift  of  integration  path  will  involve  a 

horizontal  singular  loop  integral  about  one  or  other  of  the  points  p  =  ^  M  as 

well  as  the  singular  Integral  along  the  hyperbola.  Here  we  find  that  only  residue 

terms  from  real  zeros  of  the  wave  function  (4)  make  a  real  contribution  to  the 

field,  with  the  principal  value  of  the  line  integral  vanishing  in  its  real  part. 

These  real  zeros  of  the  wave  function  provide  real  plane  wave  contributions  to 

the  reflected  shear  field,  they  are  naturally  outside  the  singular  envelope  of 

plane  waves  which  contains  the  effects  of  the  complex  poles,  and  these  head 

2  2  4 

wave  contributions  appear  only  between  the  cone  y2'*’  =  [*  +  (y  “h)  J*  >  0,  the 

2  ^ 

plane  y  =  0  and  the  tangent  planes  =  *  Mx-(y-h)(l-M  )*  . 
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Thus  the  reflected  P  field  in  y  <  0  Is 


4,  = 


{Y2T^-M^[X^  +  (y-h}^I}^ 


(5a) 


for 


+(y-h)^J^  =  Y^'^x  +  iCy 


-h,{ 


2  2  -.2,2. 
Y,t  -M  [x  +(y 


when 

YjT>(xS(y-h)V>0  . 


The  part  of  the  reflected  shear  field  inside  the  reflected  shear  cone  is 


2ir 


,  2  2  2  ,  .  ,2,^ 

[y2‘^  -(y-h)  1^ 


-  R1  < 


RlPj) 


■I' 


(5b) 


for 


P2[*^  +  (y"h)^]*  =  Y2^*+  l(y-h)[Y2T^-x^-(y-h)^]^ 


when 


Y2T>(x^+(y-h)^]*  >0 


The  associated  head  wave  contribution  is 


ill  = 


2ir 


(x^+(y  +  h)^-Y^^^ 


Im 


Ag(p)i+Bg(rti+  Cg(p)lc 
R(P) 


P  =  P 


H 


(n>^-PH)^=UPH'-n>^)^ 


(5c) 
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where  p„  is  that  value  of  the  roots 

Y2’'*  (y-h)^-  Y2‘r^]^ 

xS(y-h)^ 

2  2  — 

which  lies  between  the  point  p  =  x[x  +  (y  -  h)  ]  ^  and  the  origin. 

Now  notice  the  following  property.  In  any  given  section  of  the  transverse 

plane,  the  reflected  dilatation  cone  does  not  appear  below  the  free  surface, 

with  its  interior  field  contributions,  until  t  =  h/y^  .  Similarly,  the  reflected 

shear  cone  appears  below  the  surface  when  T>h/Y2  •  On  the  other  hand,  the 

2  i. 

planes  -  (y  -  h)  (1  -  M  )  ^ ,  which  are  tangent  to  this  reflected  shear 

cone  and  which  are  the  fronts  of  the  head  wave  contribution  can  not  appear 

2  -A 

below  the  free  surface  until  y^r  =  h{l  -  M  )  ,  when  p  =  *  M  is  a  double  zero 

of  (4)  .  There  is  therefore  a  period  when  the  reflected  shear  field  is  present 
without  a  head  wave  contribution.  The  two  stages  in  the  growth  of  this  field 
are  shown  in  figure  (2)^ 

The  remaining  terms  in  the  Integral  (1. 17)  are  more  difficult  to  evaluate. 

This  is  because  the  refracted  wave  surfaces  which  are  envelopes  of  the  plane  P 

2  2-^  2  ^ 

waves  of  form  y^T-pxf  y(M  -p  )*-h(l-p  )^  and  of  the  plane  S  waves  of 

2  -  2  2 

form  y^T-px+i^l-p  )*-h(M  -p  are  not  of  circular  section,  and  hence 
may  not  be  linked  with  simple  analytic  curves  in  the  p-plane.  What  is  clear, 
however,  is  that  the  locus  of  complex  zeros  of  either  of  these  wave  functions  may 
be  determined,  these  zeros  will  appear  in  conjugate  pairs  since  we  are  considering, 
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for  the  moment,  only  real  values  of  ,  and  therefore  the  only  contribution 
from  the  new  integration  path  comes  from  the  conjugate  zeros  of  the  wave 
function.  The  possibility  of  real  plane  wave  contributions  to  these  refracted 
fields  does  not  arise,  since  with  M<  [pi  <  1,  this  being  the  range  of  real 
values  of  p  where  we  have  seen  the  possibility  of  contributions,  neither  of  the 
refracted  wave  functions  is  real. 

This  completes  the  discussion  of  the  scattered  field  in  the  case  of  fully 
supersonic  steady  motion,  except  for  the  following  remarks.  When  h  =  0,  the 
distinction  between  the  singular  reflected  and  refracted  wave  surfaces  vanishes; 
they  are  all  of  circular  section  in  the  transverse  plane  and  all  residue  contributions 
may  be  written  down  explicitly.  Likewise  if  we  restrict  attention,  for  h  #  0,  to 
the  free  surface,  the  wave  functions  for  the  refracted  field  lose  one  of  the 
radicals  that  makes  calculations  complicated,  and  again  all  contributions  at  the 
surface  may  be  evaluated  explicitly. 

For  subsonic  values  of  the  source  velocity,  a  great  deal  of  simplicity  is 
lost.  As  the  source  velocity  a  is  reduced  from  an  infinite  value,  the  point 
p  =  M  migrates;  for  a  >  c^,  M  decreases  In  value  firom  m  *  ^  zero, 

it  becomes  negative  imaginary  when  Cj  >  a  >  c^,  and  it  becomes  real,  approaching 
the  value  1.  From  above,  when  C2  >  a-*  0  .  The  parameter  L  merely  changes 
from  real  when  a  >  Cj  to  positive  imaginary  when  a  <  Cj,  while  the  solution 
p  =  R  of  equation  (3)  has  the  complex  behaviour  of  y  as  long  as  the  source 
velocity  is  greater  than  Cj^  .  When  a  <  Cj^,  it  might  be  expected  that  like  the 


-32- 


#386 


parameter  M,  R  is  real,  but  with  R  >  1;  however,  with  the  branches  chosen, 
the  function  R{p)  of  equation  (3)  is  not  able  to  vanish  in  the  correct  range. 

This  merely  means  that  for  sources  which  move  steadily  with  a  velocity  below 
that  of  Rayleigh  waves,  the  function  R(p)  has  no  zeros,  and  thus  there  is  no 
possibility  for  Rayleigh  waves  to  be  set  up. 

Apart  from  this  point  which  is  given  special  mention,  there  are  no  surprises 
in  the  form  of  singular  surfaces.  The  singular  dilatation  surfaces  vanish  when 

^  the  singular  shear  surfaces  vanish  when  a  <  C2>  wedge-shaped 

Rayleigh  singularity  present  on  the  surface  for  a  surface  source  vanishes  if 
a  <  Cj^  .  The  other  parts  of  the  calculation  are  complicated,  both  by  the  fact 
that  the  Cauchy  principal  value  of  the  integrals  on  the  locus  of  zeros  of  the 
wave  function  is  no  longer  pure  imaginary,  and  by  the  fact  that  shifting  the 
integration  path  involves  the  calculation  of  loop  integrals  about  the  branch  points 
p  =  *  M  .  Given  these  integral  contributions,  one  might  do  just  as  well,  in 
general,  by  evaluating  the  original  integrals  on  the  formal  Integration  path.  If, 
however,  there  is  special  interest  in  fields  close  to  points  and  surfaces  of 
singularity,  then  the  shifting  of  path  is  useful,  because  the  complex  zeros  of 
the  wave  functions  appear  to  provide  as  residues  the  most  singular  part  of  the 
field  near  the  moving  source,  while  the  real  double  zeros  of  the  wave  functions 
will  give  Information  about  field  behaviour  near  the  wave  envelopes. 
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Sectlon  2;  Surface  displacements. 

Given  the  Integrals  (1. 16)  and  (1. 17)  as  general  representations  fear  the 
velocity  potentials  and  ,  we  may  calculate  both  the  velocity  and  the 
displacement  field  vectors  without  difficulty,  by  carrying  out  the  appropriate 
differentiation  and  Integration  processes  within  the  Integral  sign.  Since  the 
displacement  on  the  free  surface  Is  the  physical  quantity  of  most  Interest,  we 
shall  restrict  attention  to  the  vicinity  of  the  plane  y  =  0  . 

With  the  source  burled,  the  potentials  of  Interest  are  obtained  by  adding 
the  Integrals  {l.l6a)  and  (1.17)  .  The  velocity  vector  v  Is  V«|>+Vxi|i  ; 
this  Implies  that  the  velocity  field,  linked  In  part  with  shear  and  In  part  with 
dilatation.  Is  homogeneous  and  of  degree  -2  In  space  and  time  variables.  The 
corresponding  displacement  vector  u  Is  therefore  represented  by  homogeneous 
wave  functions  of  degree  -1  . 

For  the  burled  supersonic  source,  the  surface  displacement  Is  reduced  to 
three  torms.  These  are  given  by  the  equation 


-RlJ  2(ULV)F{2CpitLMk)(l-p^)^tI(l-LV-2p^)) 
Y,  R(P)  [y?t  -  PX  ”  -  p^)^] 


(6a) 


-00 


.R1 J  4(pC-LMA){(pL4-  LMldCl-L  V-2p^-2(M^V}^J-l(U  L^M^)(M^-p^)h 

Y2  *<P)  IY2'*’  ”  P^  ~  h(l-p^)^] 


(6b) 


2l[LMB-C(l-y)^]t2k[A(l-p^)^-Bp] 
Y7(1“P^)*[y,t  -  px  -  h(l-p^9^] 


dp 


(6c) 
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The  integral  (6a)  is  derived  from  the  P  waves  produced  by  the  source,  while 
the  other  two  integrals  are  derived  from  the  S  waves  produced  by  the  source. 

The  integral  (6b)  is  seen  to  contain  the  function  R(p)  in  the  denominator, 
while  the  integral  (6c)  does  not  contain  this  factor.  The  significance  of  this 
property,  to  be  brought  out  later  in  more  detail,  is  that  the  Integral  6c  makes 
no  Rayleigh  wave  contributions  of  any  kind,  although  some  effects  of  this  kind 
are  produced  by  (6.b)  . 

When  the  source  is  taken  on  the  surface,  the  displacement  field  has  to  be 
derived  from  the  limiting  forms,  as  h-^0,  of  the  potentials  (1. 16b)  and  (1.17). 
At  this  point  apparently  the  most  concise  way  of  defining  the  surface  displacement 
is  to  take  the  limit  as  h  -►0  of  the  integrals  (6)  and  then  to  add  the  correcting 
integrals 


oo 


2  “ 


RlJ  (PL^^LMk)(F-F)4I(M  -p)nF  4F) 

V2(Y2  (Y2T-Px)(M^-pV 


(7a) 


+  mf  LKC*+C)(l-p^)nLM(B*-B)]4l[LM(A-A*)-p(C-C*)]4kIp(B-B*)-<lV)^(A4A*)] 

YjtYz'f 

^  ^  (7b) 

this  contribution  is ,  however.  Identically  zero. 

When  the  source  velocity  a  is  taken  in  the  intermediate  range  Cj>a>c2  , 
the  velocity  parameter  M  takes  negative  Imaginary  values  as  the  transverse 
velocity  of  P  waves  becomes  imaginary.  While  this  is  the  only  formal  change 
to  be  made  in  the  integrals  (6)  and  (7),  it  will  be  found  that  there  are  line 
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integral  terms  to  complicate  the  solution  as  well  as  residue  terms.  When  the 
source  velocity  is  fully  subsonic,  besides  the  changes  in  the  velocity  parameters 
which  follow  the  change  of  both  the  transverse  velocities  from  real  to  imaginary 
values,  the  formal  integration  path  is  changed  from  the  real  to  the  imaginary 
axis,  and  an  extra  imaginary  constant  i  is  present  as  a  factor  in  the  numerator, 
just  as  in  the  integrals  (1.10b)  and  (1.12b)  . 

In  the  fully  supersonic  case,  with  F  a  real  function  of  p,  the  Integral 
(6a)  is  the  sum  of  conjugate  residues  at  complex  zeros  of  the  wave  function, 
at  any  rate  in  the  case  h  ^  0  .  The  limiting  form  of  this  sum  as  h  -*■0  must 
also  include  a  residue  contributions  from  one  of  the  Rayleigh  zeros  p  =  A  R  ^ 
because  in  this  case  one  of  these  points  lies  on  the  locus  of  zeros  of  the  wave 
function. 

Thus  for  the  buried  source,  the  surface  displacement  linked  with  the  primary 
dilatation  field  is 


-417- 


1 


2  i 

P  I 


2.  .2 


2,  2  2  2  .2,1 

c2(yiT  -X  -h  y 


R1 


(M  -p  )M2(pL+LMy(l-p  )  H1(1-L  M  -2p‘')]F 
F(P) 


(8a) 


P  =  P 


1 


2  2,2 


for 


Y-TX-iMh[Y  T  -h  -x  ] 
^  -  -  ^ 


2  u2 

X  +  h 


and  with  Ixl  <  (y^t  -h  )^  .  For  the  surface  source  we  take  the  limit  of  the 
expression  (8a)  as  h-^0,  and  the  additional  Rayleigh  residue  contribution  is 


4ir 


2.  .2 


^  j  [2itPL+LM10(p  -1)^-K1-LM  -2p  )ijF 
(dR/dp)  sgn  x 


M[t-Ix1/yj^]  ,  (8b) 


p  =  R  sgn  X 


(with  the  radicals  in  F  taking  positive  imaginary  values). 
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r  2  2  2  ,  2  i 


T  Rl< 


(pC-LMA)(l-p^)^{[l-LV-2p^-2(lV^-p^)^Kpl-HMk)-J(H-LV)(M^-p^n 


R(P) 


P=P: 


(9a) 


with 


P2  = 


V^TX  -  lh(Y^T^  - 


2  ^2 
X  +  h 


2  2  2  ^ 

When  Ixl  >  (72’’^  “  ^  >  the  real  zeros  of  the  wave  function  of  (6b)  give  the 


head  wave  contribution 


— SZsgnx  ^  ^^J(pC-LMA)(l-pV^pL-HLMjO[l-L^M^-2p^-2(I^-p^)^fl~p^y^1-Ul4-L^AK/-D^l^l 
{x^h^-Y^T^)^  I  Y2  R(P)  J 

P=Ph 

(9b) 


2222^22  1 

with  Pjj  =  [Y2TX-hsgnx(h  +x  -y^t  )2]/(x  +  h  )  and  with  the  radical  {M^-p^)2 

taking  the  specific  value,  l(Pjj-M  )*,  at  this  point.  It  is  implicit  that  the 

head  wave  contribution  only  appears  when  Ipjjl>  M;  this  limits  the  expression 

for  the  surface  displacement  to  the  range  (y^'*’^  - h^)^  <  [x]  <  YjT  -  h(l  -  /M 

In  addition  to  the  contributions  from  9a  and  9b,  taken  in  the  limit  as  h  -  0  , 

the  field  for  the  shear  source  the  surface  has  as  a  contribution  from  the  Rayleigh 

pole,  the  singular  displacement 


(pC-LMA)£pl_-fLMj^[l-L^l/-2p^  +2(p^"l^)*(p^ }  | 

(dR/dp)  sgn  x  f  1x1/ 


(9c) 
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for  p  =  R  sgn  x,  with  positive  imaginary  values  for  the  radicals  in  the  bracket 
pC  -  LMA  . 

These  representations  are  complete.  Special  notice  may  be  taken  of  the 
singular  contributions  for  the  surface  source  field  due  to  the  singular  point 
p  =  R  sgn  x;  in  general  we  find  both  a  localised  singularity  in  the  form  of  a  delta 
function  (as  in  (8b)  and  (9c),  and  a  singularity  (from  (8a)  and  (9a)  )  which 
is  locally  antisymmetric  about  the  point  jxj  =  Yj^t  .  These  contributions  are 
not  present  in  the  case  of  the  buried  source  because  the  integration  path  does 
not  then  pass  through  the  points  p  =  *  R  .  However,  when  we  evaluate  the 
residues  (8a)  and  (9a)  at  any  point  within  the  circles  |p  -  R  sgn  xj  =  R  -  1  , 
we  may  replace  the  function  R(p)  by  its  Taylor  expansion  with  the  leading  term 
of  O(p-Rsgnx)  .  For  a  given  x  and  h  the  minimum  value  of  this  factor 
occurs  when  =  Rlxj;  this  minimum  \elue  is  small  if  h<<  x,  and  it  is 
found  that  a  surface  disturbance  whose  amplitude  is  0(x/h}  will  pass  a  fixed 
point  with  the  velocity  of  Rayleigh  waves,  with  its  peak  seen  for  =  [xj  , 

Without  numerical  work  it  is  not  possible  to  say  whether  this  maximum  is 
noticeable  when  h/x  is  not  small.  It  is  however  implied  in  the  formulae  (8a) 
and  (9a)  that  the  quantities  and  p^  are  complex;  the  application  of  Taylor 
series  to  the  residue  (8a)  is  therefore  restricted  to  the  range  h<  |xI(R  -  1)^  . 
These  limits  have  been  mentioned  elsewhere  (Ewing  et  al. ,  1957),  and  they 
differ  from  the  empirical  resvilt  given  by  Pekeris  (1957). 
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Part  3;  The  transient  moving  source  In  an  elastic  solid  with  a  free  plane  surface. 


Section  1:  The  simplification  of  double  Integral  solutions. 


For  the  steady  supersonic  source  we  have  given  a  fairly  detailed  description 

of  the  process  by  which  Integrals  of  singular  plane  waves  may  be  reduced  to 

simple  residue  calculations.  For  the  double  integrals  associated  with  a 

transient  source,  the  evaluation  is  much  more  difficult;  with  two  integration 

parameters  there  is  no  unique  choice  of  complex  paths,  for  each  parameter,  on 

which  the  wave  function  takes  real  values. 

Fortunately,  while  we  are  dealing  with  propagation  in  an  isotropic  medium, 

we  are  able  to  determine  a  large  amount  of  information  in  the  manner  described 

below.  The  simplest  situation  involves  the  reduction  of  the  point  source  fields 

in  the  absence  of  a  boundary.  In  the  integrals  (1.  21b)  or  (1.23)  we  have 

specific  integrals  for  the  dilatation  potential  of  a  fixed  transient  source.  In 

“l  “3 

these  integrals  we  are  free  to  shift  the  'slowness'  parameters  (  — ,~,  — ) 


around  the  surface  of  a  sphere  of  radius 
2  2  -  2  “ 

and  a^=(m  -p)^(l-q)^,  so  that  a 


m  =  Cj/c^;  we  may  put  aj=(m^-p^}*q 
=  p,  and  we  are  to  evaluate  the  integral 


Rijjr  ^-T - 2  ^ 2^i *  ^ — 2^ - 

)^{c2t-(m  -p  )^[qx+(l-q  )^z]-pysgny} 
where,  in  the  case  of  a  source  of  unit  strength, 

F  =  l/8Tr^  , 


(2a) 


#386 


-39- 


and,  in  the  case  of  a  point  force 

F  =  -  {(m^  -  [qX  +  (l -q^)^Z]+ pY  sgn  y  }/8ir^C2  .  (2b) 

For  X  =  r  cos  0,  z  =  r  sin  0,  we  may  choose  a  complex  path  on  which 

2  i 

q  =  cosh(w  +  i0)  sgn  z,  (1-q  =  -i  sinh(w+i0)  sgn  z  , 

so  that 

2  i 

qx  +  (1  -q  )*z  =  r  sgn  z  cosh  w  . 

On  this  path  the  integral  (1)  has  the  form 

4,  -  -  /dw  R1  / cosh (w+ 10)  sgn  z]  +  F[p,  cosh(w-i0)  sgn  z] 

-00  -00  [c^t  -  (m^  -  p^)^r  cosh  w  sgn z  -  py  sgn  y 

As  long  as  the  function  F(p,  q)  is  a  real  function  on  the  real  axis  of  both  variables 
when  Iql<  1  and  Ipi  <  m,  (as  it  is  in  the  specific  cases  given  in  equation  (2)), 
there  is  a  real  contribution  to  the  integral  (3)  only  from  residues  at  conjugate 
poles  of  the  integrand.  The  integration  path  for  p  is  taken,  just  as  in  the 
simpler  case  of  steady  motion,  along  the  locus  of  complex  zeros  of  the  demonlmator 
for  fixed  values  of  x,  y,  z  as  t  is  varied.  The  principal  value  of  the  integral 
along  this  path  is  easily  shown  to  be  the  difference  of  conjugates  and  hence 
imaginary,  but  the  residue  contributions  from  the  zeros  of  the  denominator  give 
a  real  expression.  Thus  then  the  complex  zero  in  the  fourth  quadrant  of  the 
p-plane  is  at  the  point  p  =  Pj  where 
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^1  = 


2  2  2  2  2 

c^ty  sgn  y-lrm  cosh  w  (c^t  -y'-r  cosh  w) 

2  2  2 
y  +  r  cosh  w 


(4) 


z  z  - 

this  being  a  possibility  only  if,  for  c^t  >  R  =  (r  +  y  )^  >  0,  w  is  restricted 

2  2  2  - 

to  the  range  for  which  cosh  w  <  (c^t  -y  )Vr  .  The  evaluation  of  residues  at 
Pj  and  its  conjugate  point  gives  the  result  that 


?  O  1 

[(m  -p  )^{F[p  ,  cosh(w+ie)sgn2]  +  F[p  , 

<j)  =  -  2ir  J  2  2  2  2  2 

0  nfi(c^t  -y  -r  cosh  w)^ 


cosh(w  - 10}  sgn  z]}] 
-  dw 

(5) 


2  2  Z  — 

with  Wj  =  arccosh(Cjt  -  y  )^/r  .  This  is  the  only  contribution  for  the  fixed 
transient  source.  With  F  given  as  in  (2),  the  finite  integration  can  be  carried 
out  explicitly.  Note  that  the  only  contribution  in  these  integrals  is  confined  to 
the  interior  of  the  wave  surface  R  =  c^t,  which  is  a  singular  surface  fcM:  the 
integral  (1)  because  it  is  the  envelope  of  the  singular  plane  waves  given  in 
that  expression.  The  same  treatment  leads  to  the  potential  for  a  source  moving 
with  subsonic  velocity,  the  whole  of  the  field  being  contained  within  the 
spherical  wave  surface  R  =  c^t  .  Thus  if  we  take  the  case  of  a  transient  source 
which  moves  along  the  y-axis  with  a  velocity  a  we  have  to  examine  the  integral 


00 


i  c^  F(Pjq) 


2-  22-  2- 
-oo(i-q  )2[c2-asgnyp]{c2t-(m  -p  )^[qx  +  (l-q  )®z]-py  sgny} 


dp  dq 


(6) 
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When  a  <  the  pole  P  =  sgn  y/a  can  be  identified  with  a  specific  value 
of  Pj  only  when  r  =  0  and  y  =  at,  and  this  merely  gives  us  a  singular  value 
at  the  source  point  for  the  complete  field  in  the  form 


w 


Zir/ 


1 


c^  dw 


z  ,  2  2  2  2  .2  ,4 

0  [c^t  -y  -r  cosh  wj^ 


Rl< 


/  2  2  i 

(m-p^ 


(c^-a  sgnyp 


ry  P^^l*  cosh(w+10)sgn z]  +  F[Pj,  cosh(w-10)sgn  z] 


(7) 


When  a>Cp  the  same  pole  lies  between  the  branch  points  p=*m,  and  it 
coincides  with  a  real  zero  P  =  P^  of  the  denominator  of  the  Integrand,  where 


P2  = 


2  2  2  ..  w 

Cjtysgny  *  mrcoshw(y  +  r  cosh  w  -  c“t“) 


2.2,4 


2  2  2 
y  +  r  cosh  w 


The  point  about  these  real  zeros  is  that  neither  is  capable  of  generating  a  real 
residue  contribution  to  the  integral  (3)  or  (5)  unless  the  remaining  factc^  of 
the  integrand  is  non-real,  or  singular,  at  the  point  concerned.  There  is  an 
Important  choice  to  be  made  here.  For  the  integral  (6)  only  the  smallor  value  of 
can  make  a  contribution  in  this  manner;  this  is  the  value  which  lies  to  the 
left  of  the  locus  of  complex  zeros  for  a  given  y,  r.  Thus  when,  for  y  >  0  and 
R  >  Cjt,  we  are  able  to  find  a  p^  with  the  value  c^/a  and  we  may  take  the 
residue  at  this  point  in  the  negative  sense  to  obtain  the  additional  term 


2"Y,  « 


*m00 


(1-q^)^  (Yj  (t  -  y/a)  -  qx  -  z(l  -  q*")  *] 


.Uc/SL  qL 


—  me  —  yi  1  —  a^\  *  1 


dq 


(8) 
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This  term  is  clearly  a  steadily  moving  supersonic  field  trailing  behind  the  source 
point  and  contained  within  the  singular  conical  front  defined  by  the  equation 

this  cone  being  a  tangential  surface  to  the  sphere  R  =  c^t  . 

The  reason  for  choosing  the  pole  to  lie  to  the  left  of  the  complex  integration 
path  is  that  this  restricts  the  conical  ccntribution  to  the  region  between  the  wave 
surface  R  =  c^  t  and  the  source  point.  If  the  other  real  zero  of  the  wave  function 
is  involved,  it  produces  the  complementary  conical  field  trailing  behind  the 
sphere  R  =  c^  t,  but  only  for  y>  0  .  This  is  not  permissible  because  it  involves 
discontinuity  across  the  fixed  surface  y  =  0,  and  this  plane  is  not  a 
characteristic  surface  of  the  system. 

To  make  the  same  point  more  explicit,  we  should  state  that  in  the  integral 
(1),  the  formal  integration  path  for  q  is  along  the  real  axis  except  for  deformations 
below  q  =  -1  and  above  q  =  1,  and  that  for  p  is  also  along  the  real  axis,  but 
passing  below  the  point  p  =  -  m  and  above  p  =  m,  above  the  smaller  value  of  p^ 
and  below  the  larger  value  of  p^  . 

The  importance  of  the  choice  of  this  path  is  made  more  evident  if  we  change 

2  2*^  2  2  4 

variables  of  integration,  putting  P  =  (m  -p  )*  and  (m  -P  )*  =  p  .  Then  with 

_  ^  i  P  P[(m^-P^)S  q]  c^  dp  dq 

^  “  R1 24224  2  2  i  2  i  ?74. 

-00  (1-q  )»(m  -P  )*[c2-a(m  -P  )*sgny][c2t-P(qx+(l-q  J*z)-(m^-P^)*  y  sgny] 


(9) 
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the  supersonic  contribution  of  this  integral  is  found  to  be  linked  with  the  real 
zero  of  the  wave  function  which  lies  to  the  right  of  the  complex  Integration  path 
and  which  coincides  with  the  point  P  =  c^/yj  0  •  For  the  purpose  of 

deciding  the  sense  of  integration  at  complex  poles  we  may  take  the  formal 
integration  path  for  P  along  the  real  axis,  below  the  point  P  =  -  m,  above 
the  point  P  =  m,  and  below  the  smaller  and  above  the  larger  of  the  real  zeros 


P  = 


2  2  2  2  2 
c^tr  cosh  w*  my(y  +r  cosh  w-c^t  ) 


2  2  2 
y  +  r  cosh  w 


The  next  point  of  difficulty  in  the  evaluation  of  the  scattered  field  is 
connected  with  the  presence  of  the  characteristic  Rayleigh  function 

R  =  (1  -2«^^  -  +  4  +  a^) 

which,  because  it  has  its  own  zeros,  can  contribute  its  own  singularities  to 
parts  of  the  field.  If  we  choose 

2  i 

=  pq,  #3  =  p(l-q 

with 

a=(m^-p^)2,  p^(l-p^)2 

then 

R  =  R(p)  =  (1  -  2p^)^  +  4p^(l  -  p^)^(m^  “  P^)*  , 

and  this  function  is  easily  recognised  to  have  zeros  at  the  points 

P  =  *  °2/°R  » 

c  1^  being  the  velocity  of  Rayleigh  waves. 
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The  characteristic  function  R  only  appears  in  the  expressions  1. 27  for 
the  scattered  field  in  idie  region  y  <  0  .  The  complex  integration  path  for 
these  expressions  will  only  pass  through  these  Rayleigh  zeros  when  foe  a  surface 
source  with  h  =  0  the  point  of  observation  is  also  on  the  surface  with  y  =  0  . 
Only  in  this  case  do  we  find  a  singular  surface  wave  associated  with  the 
transient  source.  For  the  source  which  is  initiated  at  a  depth  h,  there  will 
be  found,  when  h/r  <<  1,  a  disturbance,  which  travels  with  velocity  Cj^  , 
whose  amplitude  is  0(r/h)  .  This  point  arises  just  as  in  the  case  of  steady 
source  motion.  Even  when  the  buried  source  moves  towards  the  surface  after 
it  is  set  up,  there  is  no  singular  surface  wave  until  the  source  arrives  at  the 
surface,  and  then  only  because  whatever  may  happen  to  the  source  there  is 
bound  to  be  a  transient  process  at  this  moment  which  will  set  up  its  own  spherical 
wave  surfaces. 

We  have  been  discussing  the  various  difficulties  which  arise  in  the 
evaluation  of  the  dilatation  potentials.  The  same  difficulties  also  arise  in  the 
discussion  of  the  shear  potentials,  but  there  is  one  more  special  situation, 
linked  with  the  presence  of  head  waves. 

Head  waves  arise  in  the  reflected  shear  field  part  of  the  integral  (1. 27b}  . 

If  we  make  the  choice  of  and  as  given  in  the  equations  (10),  it  is 
clear  that  although  the  wave  function  contains  only  the  branch  points  p  =  *  1  , 
the  remaining  factors  of  the  integral  (specifically  the  Rayleigh  function  R  ) 
contain  the  branch  points  p  =  *  1  and  p  =  ^  m  .  We  have  now  the  possibility 
that  real  zeros  of  the  wave  function  can  appear  in  the  range  m<  |p(  <  1  ,  where 
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the  Rayleigh  function  is  complex,  and  from  these  zeros  we  can  find  real 
contributions  to  the  potential.  The  head  wave  contribution  arises  because  the 
refracted  P  field  travels  faster  along  the  surface  than  the  associated  shear 
field,  and  is  incapable  of  satisfying  the  boundary  conditions  there  by  itself. 

The  support  for  the  head  wave  field  is  thus  on  the  plane  y  =  0,  and  therefore 
it  is  real  zero  of  the  wave  function  which  lie  between  the  branch  point  p  =  m  sgn  z 
and  the  complex  locus 


P  = 


2  2  2  2  2 

c^tr  cosh  w  sgn  z  ^  iy  (c^t  -y  -r  cosh  w) 


2  2  2 
y  +  r  cosh  w 


that  is  the  zero 


%  ■ 


22  2  22- 
[c^trcoshw  +  y  (y  +r  cosh  w  -  c^t  )^J 

2  2  T 

y  +  r  cosh  w 


sgn  z 


which  gives  a  real  residue  contribution,  with  the  range  of  w  restricted  to  keep 
2  2 

1p„1  >  m  when  (y  +  r  )^  =  R  >  c-t  . 

ri  L 

There  remain  only  computational  difficulties  in  the  evaluation  of  the  total 
field  due  to  a  transient  source.  The  primary  source  fields  and  the  reflected 
fields  can  all  be  written  down  explicitly  as  finite  integrals  with  respect  to  w  , 
the  ease  with  which  these  integrals  can  be  found  being  linked  with  the  simple 
nature  of  the  associated  singular  wave  fronts  (they  are  either  spherical  or  conical). 
The  refracted  wave  fronts  are  not  so  simple  to  define  without  the  use  of  parameters, 
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but  because  the  complex  zeros  of  the  wave  functions  arise  in  conjugate  pairs, 
there  is  no  doubt  about  the  reduction  of  the  p-integration  to  a  residue  calculation; 
these  zeros  however  are  found  as  solutions  of  a  quartic  algebraic  equation. 

Finally  it  will  be  repeated  here  that  the  structure  of  the  scattered  field  is 
determined  by  the  singular  surfaces  for  the  Integrals  (1.27),  and  these  surfaces 
are  determined  by  the  condition  that  the  denominator  and  its  derivatives  with 
respect  to  both  Integration  parameters  vanish  together.  There  is  no  difficulty 
in  this  calculation. 
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Section  2:  Calculations  of  surface  displacements  for  a  vertical  force. 

We  may  derive  expressions  for  the  surface  displacement  due  to  the  general 
point  source  with  a  homogeneous  potential  of  degree  -1  .  The  potentials 
(1.  26)  and  (1.  27)  are  first  differentiated  to  form  the  velocity  and  then  integrated 
with  respect  to  time  to  form  integrals  for  the  surface  displacement.  For  the 
transient  source  which  appears  at  a  depth  h,  we  find  three  distinct  integrals; 
these  are,  first 


00 


i  rr  J2a[a^  i_+ Q-^jy +  (l-2«^  -  2a^)J.} 

Up  =  jj  - - rr—r - :: - r— n - da. 


2  -00 


[c^t-o^x-a^z-ah] 


dflf. 


(lla) 


second 


^SR  " 


4  ^  -Aa  ]{{o  i_+ar  kKl-2a^-2a^-2flfp)-la} 

^2  R(«i,«3)[c2t-ajXa2Z-ph] 


1  *^"3  ’ 


(lib) 


and  last 


2  _  “  1  V^([Ba3  -  CP]  1. 1  [tf  - 


-00 


P[c2t  -  a^x  -a^z  -  ph] 


d«i  d«3 


(lie) 


Of  these  integrals,  the  first  is  linked  directly  with  the  primary  P  field;  the 
others,  linked  directly  with  the  primary  S  field,  differ  in  that  the  last  is  a 
simple  horizontal  displacement  field  which  contains  neither  head  wave  nor  surface 
wave  contributions,  while  (11b)  will  be  shown  to  contain  both  these  contributions 
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when  appropriate.  (The  reason  for  the  simplicity  of  (11c)  is  that  the  potentials 

2  2  2 

A,  B  and  C  may  contain  neither  the  radical  a=  (m  “*3^* 
characteristic  Rayleigh  function  R(aj,  “3)*)  These  integrals,  with  h  =  0,  are 
correct  expressions  for  the  surface  displacement  when  the  transient  source  is  set 
up  on  the  surface. 

The  pure  dilatation  source,  with  F  constant,  A  =  B  =  C  =  0,  is  the 
simplest  case  to  evaluate,  but  since  a  head  wave  disturbance  is  not  present 
for  the  buried  source,  certain  aspects  of  the  process  of  evaluation  of  the 
Integrals  (11)  are  not  brought  out.  The  simplest  case  we  can  take  for  a  detailed 
examination  is  that  of  the  transient  vertical  force,  acting  at  a  fixed  point.  For 
the  vertical  force  Y  we  have  the  potentials 


F  =  -(m^-p^)^  Y/Bir^c^,  B  =  0  and  Ai.+ Ck  = -Yp[(l-q^)^i.-qkJ/8ir^C2  , 

(12) 

2  ^ 

with  Vj  =  V2  =  1,  these  being  the  forms  when  cj  =  pq,  ofj  =  p(l-q  )*  . 

The  integrals  to  be  evaluated  are 

-P  -00  (l-q^)2R(p){c2t-p[qx  +  (l-q  )^2]-(m  -p)*h} 

(13a) 

-«>  (l-q^)*R(p){c2t-P[<JX  +  (l-q^)*  z]-(l-p^)*  h} 


,32 
2ir  c^ 


(13b) 
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and 


oo 


>.  3  2 

4.  Cj 


2  i  2  ^  2  i 

-*  (i-q  )^(o,t-ptqx+ (i-q  )‘2)-(i-p  )“h} 


■ujr 


dp  dq 


(13c) 


Now  take  the  integration  path  for  q  into  the  complex  plane;  with 

2  ^ 

x=rcos0,  z=rsin0,  q  =  cosh(w  +  i0)  and  (1 -q  }^  = -1  sinh(w  +  10)  ; 

(with  a  restriction  to  z  >  0  in  order  to  simplify  the  writing  down  without  losing 
any  generality)  we  find  that  on  this  path  the  contribution  from  w  <  0  is  the 
complex  conjugate  of  the  contribution  for  w  >  0  .  Hence  the  q-integratlon  may 
be  reduced,  in  part,  to  an  integration  for  positive  w  ,  with  the  contributions 

„(1)  =  -  ^  r“ dw  R1  /  ani-pV y-lL  dp  ,  (14a) 

— ^  0  -»  R(p)  [c2t-prcoshw-(m  -p)*h] 


where  r 


u(i)^  ^  j^dw  Ri  /  -  2(1  - p^)^ - p^fllrC^^  - p¥i> 

R(p)[c2t-prcoshw-(l-p^)2  h] 


SR  3  2  „ 

V  C2  0 


-00 


{14b) 


and 


*  ^  p  cosh  w  r 

/  dw  Rl  J  - *- 


2  - 

-00  c2t-prcoshw-(l-p  )*  h 


dp 


following  simply  from  the  Integrals  (J3)  . 


(14c) 
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The  integration  with  respect  to  p  may  be  carried  out  in  the  usual  manner. 

Thus  in  (14a),  the  only  real  contribution  to  the  Integral  comes  from  residues  at 

the  conjugate  complex  zeros  of  the  wave  function,  and  these  complex  zeros 

2  2  Z  - 

are  only  present  under  the  restriction  that  coshw<(Cjt  -h  )*/r  .  Thus 


.CD 


arccosh  (c^ t^  -  h^]^/r 


2  ,2  ,2 
h  -r  cosh  w)' 


T  Rl{G^(Pj)} 


where 


^1  = 


c^  tr  cosh  w  -  imh  (c ^  t^  -  h^  -  r^  cosh^  w) 

2  2  2 
h  +  r  cosh  w 


2  X  2  2  ~ 

where  both  the  radicals  (1-p  and  (m  ~p  have  positive  real  parts  when 


p  =  Pp  and  where 


G  ini  =  -  P^)  11(1  -  2p^)  2r.  P(1  -  P^)^  cosh  w]  Y 

m  R(p) 


2  2  2  - 

This  Integral  is  only  present  when  O-^r  "^(c^  t  -h  .  To  reduce  it  to  a 
simpler  form,  we  may  make  the  real  transformation 


2  2  .2  2  .2  ,  2  2  .2  2,  ,  2^ 

Cjt  -h  -r  cosh  w=(Cjt  -h  -r  )  sin  \  , 


and  we  have  the  result  that 
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tt/2 

■  "  J. 


{Gj(pp} 


(15a) 


2  2  Z  ^ 

within  the  singular  P  circle  r  =(Cjt  -h  )*  . 


With 


2  2  2  2  Z  ^ 

2  c^tr  cosh  w  -  ih  (c^t  -  h  -  r  cosh  w)^ 

p  _  2  2  2 

h  +  r  cosh  w 


with  positive  real  parts  specified  for  the  radicals  and  under  the  transformation 


2*2  ,2  2.2  ,  2  2  2,  ,  2 

c^t  -  h  -  r  cosh  w  =  (c^t  -  h  -  r  )  sin  p  , 


the  shear  field  contributions  from  complex  zeros  of  the  wave  function  are 


ir/2 

f  - 

SR  22  /  C 

■  ^  ^2^  I 


(15b) 


with 


^  Y{2p^(l-p^)2coshw[l-2p^-2(l-pV(m^-p^)*£-2(l-p^)2(m^-p^)Vl} 

C^lp)  —  __ 


dnd 
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""S"  2 
—  ir 


Tr/2 

^  f- 

vl 


col^ 


(15c) 


with 

G^Cp)  =  Yp(l  -  cosh  w  £  , 

2  2  2  “ 

inside  the  singular  S  circle  r  =  (C2t  -h  . 

For  the  buried  source  with  h^O,  the  integrals  (15a)  and  (15c)  are  the 
only  contribution  to  the  surface  displacement  to  be  derived  from  the  integrals 
(14a)  and  (14c)  .  The  integral  (14b),  on  the  other  hand,  contains  both  the 
contribution  (15b)  linked  with  complex  poles  of  the  Integrand  and  a  head  wave 
contribution  as  well.  There  is  a  real  residue  contribution  from  the  real  pole 
p  =  Pjj  where 

2  2  2  2  2*^ 
c^tr  cosh  w  -  h(r  cosh  w+h  -c^t)^ 

Pjj  =  2^  2  72  ' 

h  +  r  cosh  w 

this  contribution  appears  only  when  p.,  >  m  .  This  condition  imposes  the 

2  ~ 

restriction  that  r  cosh  w  <  c^t  -  h(l  -  m  )^/m  .  With  the  change  of  integration 
variable  given  by 

r^  cosh^  w  +  h^  -  c^t^  =  (r^  +  h^  -  c^t^)  cosh^  A 
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we  find  that  for  -  h^)^<  r  <  c^t  -  h(l-m^)Vni  , 


“1 

f  ■ 

SR  2  2  < 

-  n  c^r  J 


Im  {G-CPu)} 
coshw  *■ 


(15d) 


with 


=  arccosh 


'  2i  ^ 

-  m  -h 

mCr^  +  h^-c^)"  ^ 


This  is  the  head  wave  contribution  for  the  buried  source. 

2  Z  2 

Notice  that  for  points  on  the  compressional  wave  front,  =  in(l  -h  /c^t  ) 
is  constant  and  the  integral  (15a}  takes  the  simple  value 


3 


=  Gj(Pj)/2v  c^r 


2  2  ^ 

with  w  =  0  .  On  the  shear  front  p^  =  (1-h  /c^t  is  constant,  and  the 
integrals  (15b)  and  (15c)  similarly  that  the  simple  values 


4r  = 


Ug  =  GjCPjj/Z-r  o  r  , 


and 


Him 


-54- 


#386 


with  w  =  0  .  These  simple  expressions  are  the  initial  displacements  on  the 

surface  linked  with  the  arrival  of  P  and  S  waves;  they  agree  with  the  results 

of  Pekeris  (1957)  for  the  horizontal  components  of  displacement,  but  for  the 

vertical  component  they  do  not.  (His  initial  P  displacement  lacks  the  factor 
2  2  2  - 

(l-h  /Cjt  ,  and  his  initial  S  displacement  lacks  the  factor 

Ic^  .) 


The  Integral  for  the  head  wave  is  also  simplified  at  the  front  of  the  S  wave. 

The  integrand  (15d)  reduces  to  a  constant,  but  the  integral  becomes 

logarithmically  Infinite  together  with  the  upper  limit,  A^,  of  integration. 

Pekeris  also  notes  this  logarithmic  singularity;  we  differ  in  that  he  places  it 

just  inside  the  S  wave,  instead  of  just  outside,  and  in  that  his  formula  for 

2  2 

this  logarithmic  singularity  in  the  vertical  component  lacks  the  factor  (1  -h  /c^t  ) 
For  the  surface  source,  the  integrals  (15)  have  to  be  evaluated  with  h  =  0  . 
This  is  not  a  uniform  limit  of  the  case  for  the  buried  source  because  the 
integration  path  Involves  real  values  of  p  to  the  right  of  p  =  m,  thus  passing 
through  the  Rayleigh  pole  p  =  R  .  Thus  with  the  total  surface  displacement 
field  (14)  written  in  the  more  concise  form 


<»  «_/  2  2.  i  ...  2  ,,,  2.i  2  2.i  2- 

Ri  r  H w  r  -2(l-p  )^(m  -p  )^]p  ^coshw 

^  io  R(p)[c2t-pr  cosh  w]  P 


u  =  - 


o  3  2 

2ir  c^ 


(16) 


and  with  the  p-integratlon  path  shifted  to  the  horizontal  loop  to  the  right  of  p  =  m 
we  find  a  number  of  distinct  contributions.  These  are 
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2y  [  (1  -  2p^)(p^-m^)^p[(i-2p^)J.  +  (1  -  £c^t/r] 

J  [(l-Zp^)** +16p^(l-p^)(p^-m^)](C2t^-pVj^ 

m 


when  m  <  c-t/r  <  1,  and 


u=-^  { 
-  2  2/ 

"  °2  J 

m 


^  pC^  -  (P^  -  m^)ni  -  p^)^i  c^t/r 


r/i  o  2.4^,,  4,,  2„  2  2.,,  2^2  2  2.4 

[(l-2p  )  +16p  (1-p  }(p  -m  )](c2t  -p  r  )* 


dp  (17b) 


Cjt/r 


2  2 

2Y  r  _ pCp  -m  )^1 _ 

2  2/  ,2,2  ,2,2  „4,  2  2,4,,  2  2  2  2,4 

^  ^  ^P  “1)  (P  ”«»  )  1(^2^  "P  ^  J 


dp  (17c} 


_2Y_ 

2  ,  2  2  „2  2,4 

1^02  D(R)  (02^  -R  r  )* 


(17d) 


where  D(R)  =  [dR(p)/dp]  = -8R[6R^(l-in^)+2R^(2m^-3)  +  l](l-2R^)"^  , 

p  =  R 

when  c^li>l  for  (17b),  (17c),  and  for  (17d)  . 

The  expressions  (17a}  and  (17b)  come  firom  the  contribution  to  (16)  of 
the  range  m  <  p  <  1,  (17c)  comes  from  the  range  P  >  1,  and  (17d}  comes  from 
the  single  point  p  =  R  .  No  other  real  terms  are  present.  The  Integrals  (17a} 
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and  (17c)  contain  a  singularity  at  p  =  R  if  C2t/r>R  . 

The  integrals  are  easily  Identified  with  those  given  by  Pekeris  (1955a}:  we 

differ  in  the  value  of  the  coefficient  of  the  final  Rayleigh  wave  term.  It  is  noted 

that  the  integrals  (17)  may  be  simplified;  the  vertical  component  may  be 

evaluated  explicitly  and  the  horizontal  component  reduces  to  forms  involving 

incomplete  elliptic  integrals.  This  work  has  been  done  by  Pekeris  (i955a). 

[What  is  most  interesting  is  that  the  surface  displacement  field  may  be  found  in 

its  simplest  form  by  an  alternative  integration  process  (for  (13))  involving  first 

the  calculation  of  residue  contributions  from  the  zeros  of  R(p}  due  to  ^  the 

2  “  2  2*^ 

branches  of  the  radicals  (1-p  )*  and  (m  -p  ).*,  and  then  the  evaluation  of 
the  q-integral  in  the  usual  manner.  For  the  vertical  component  the  q-integration 
is  simply  a  matter  of  calculating  residues,  and  this  is  a  simple  way  of  noting 
the  existence  of  a  lacuna,  already  seen  by  Pekeris.  The  horizontal  component, 
however.  Involves  quadratures.  The  alternative  method  of  integration  which  is 
centred  on  the  existence  of  multiple  sheets  for  the  Rayleigh  slowness  surface 
will  be  discussed  elsewhere.  ] 

Notice  finally  that  the  appearence  of  Rayleigh  singularities  far  the  surface 
source  is  accompanied  by  the  simultaneous  reduction  in  the  order  of  singularity 
linked  with  the  arrival  of  P  and  S  waves.  Where  for  the  buried  source  we 
expect  a  step  function  singularity  for  the  P  wave,  or  a  step  function  coupled  with 
a  conjugate  logaritmio  singularity  for  the  S  wave,  the  surface  source  field  has 
the  same  singularities  in  its  first  derivatives,  and  is  therefore  continuous  in  its 
displacement  components. 
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Sectlon  3;  Surface  displacements  for  fixed  sources  with  an  unsymmetrical 
displacement  field;  the  horizontal  force  and  the  couple  with  arbitrary 
orientation. 

The  value  of  the  present  method  is  that  the  evaluation  of  the  surface 
displacement  field  from  the  integrals  (11)  is  not  made  any  more  difficult  if 
there  is  no  axis  of  symmetry.  Having  just  considered  the  case  of  the  vertical 
force,  with  its  obvious  axis  of  symmetry,  we  shall  now  examine  the  case  where 
a  horizontal  force  Xi.  is  applied  at  a  point  below  the  surface  of  an  elastic  solid 
With  the  potentials 

F  = -Xpq/8  ir^C2,  A  =  0,  B  =  Xp(l  - q^)V8ir^C2,  C  = -X(l  -  p^)  V8ir^C2 

(18) 


the  integrals  11  take  the  form 


~  3  7  7  7  i  7  — 

-  c  -00  R(p)(l-q  )2{c  t-p[qx(l-q  )2z]-h(m  -p**)*} 


*  ip^q{2p(l-p^)^[cU, +(l-q^)*l^  +(l-2p^)i} 


?  i  dp  dq  , 
(19a) 


“SR" 


Zir  c^ 


*  lp^q(l-p^)^{p  [qL+  (1  -  q^)  ^  k  J[l-2p^-2(l-p^)^(m^-p^)^]  -jj m^-p^)^} 

-00  R(p)  (l-q^)Mc2t-p[qx+(l-q^)2z]-h(l-p^)2} 

(19b) 


^If 


and 


u„  = 


X 


.,3  2 
4tr  C2 


ip^q  [gLl  ^(i-q^)^  kj 


(1  -  p^)  ^(l-q^)  ^  {c^t  -  p[qx + (1  -  q^)  2 


z]-h(l-pV} 


dp  dq 


(19c) 


dqdq, 
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It  follows,  with  the  abbreviations  cosh  w  =  c,  sinh  w  =  s,  f  cos  0  +  k  sin  0 
and  £  =  k  cos  0  sin  0,  that 


Z  2  Z  ^  f  Z  2*^2  2  I 

G^(p)  m  R(p)  = -Xp  (m  -  p  |(1  -  2p  )  c  cos  0  j.+ Zp^l  -  p  )*  [c  fcos  S+ s  0sin0]|, 

(20a) 

G,  (p)  RCp)  =  2Xp^(l  -  p^)  (m^-p^)  ^  c  cos  01  +  p  [1  -  2p^  -  2(1  -  p^)  ^  (m^  -  p^)  *]  • 


2  «  2  - 
[c  xcos  0  +  s  Q  sin  0] 


(20b) 


and 


G,(p)  =  Xp  [r  cos  0  (1-p^c^)  -  0  sin  0  (1  +  p^s^)] 


(20c} 


these  being  the  three  expressions  needed  in  the  integrals  (15)  to  define  the 
surface  displacement  for  the  burled  source. 

For  the  surface  source  we  have  the  additional  term  associated  with  the  residue 
at  the  point  p  =  R  .  The  complete  contribution  to  the  surface  displacement  of 
this  residue  is  restricted  to  the  region  r  <  c_t ,  with  the  form 

_  Xt[l-2R^  +  2(R^-l)^(R^-m^)^Il 
irc2rD(R)  (c^t^-r^)^ 

The  field  due  to  the  sudden  application  of  a  couple  is  of  interest  to 
seismologists:  this  is  the  point  source  model  of  symmetrical  shear  in  the  focus 
of  tectonic  earthquakes,  with  symmetrical  displacement  in  relation  to  a  plane  of 
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rupture,  and  with  movement  of  the  sides  of  this  plane  in  opposite  directions 
(Kellis-Borok  I960}.  The  couple  is  therefore  to  be  taken  as  the  limit  of  a  double 
force  with  moment,  with  the  force  acting  in  the  plane  of  rupture. 

We  specify  that  the  couple  is  of  moment  M,  that  it  acts  in  the  plane  z  s  0  , 
and  that  the  component  forces  in  this  plane  act  at  an  angle  c  to  the  plane  y  =  0  . 
For  the  couple  of  step  function  time  dependence,  the  velocity  potential  and  the 
displacement  field  are  both  functions  of  degree  -2  in  space  and  time  variables. 

To  keep  within  the  structure  of  the  calculations  of  this  paper,  we  shall  be 
considering  the  surface  displacements  given  in  the  equations  (11)  and  these 
must  subsequently  be  differentiated  with  respect  to  t  .  It  might  be  suggested 
that  the  result  for  a  couple  is  easily  found  by  taking  the  results  for  a  force  and 
differentiating  in  the  appropriate  direction;  this  process  is  however  available 
only  when  we  have  explicit  results  for  the  field  both  on  and  off  the  free  surface. 

Thus  we  use  the  potentials 

F  =  — ~^(pqcos€  +  (m^-p^}^slne][pqsln«  -(m^-p^)^cos€]  , 


=  — 7-7:[pqsin«  -(1 -p^)^cos€] 


p(l-q  )*  Sint 

-pO-q^)^slne  • 

2  4 

(1-p  )*cose  -pqsinc 


,  (21) 


in  the  integrals  (11),  and  we  derive  the  quantities 
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Gj(p)2mc2R(p) 


p(m^-p^)*  M 


[2p(l-p^}^ci+(l-2p^)l] 


2  2  2  2  2  2  2 
sln2€[p  c  cos  20 -p  s  sin  0+p  -m  ] 

2  2  ~ 

-2(m  -p  )*cos  0COS  2« 


2  2^2-^  22*^ 

2p(l-p)*s  0sin0[2(m  -  p  )*  cos  2€  -  pc  cos  0  sin  2£] 


g^Cp)  c^RCp) 

2  Mp^(l-p^)^ 


{pci[l  -  2p^  -  2(1  -  p^)*  (m^  -  pV]  -l(m^  -  p^)*  }  • 

,,  2,4  2  2^  , 

p(l-p)'s  sin  0 Sint  cose 

+ 

2  -  2  — 

[(1-p  )^cos€  -pccos0sin£][c(l-p  }‘cos0cosc  -psint] 


-  ps^  sin  0  0  [1  -  2p^  -  2(1  -  p^)*(in^  -  p^)*]  • 

2  4 

psint  [c(l-p  }‘cos0cost  -psint]  + 

i  2  4  2  ^ 

+  s(l-p  )*cost  [(1-p  )*cost  -pccos0sint] 


and 


°3(Pl  °2 

Mp 


[pc  cos  0  sin  t 


-(1-p^)*  cos  t] 


A  2  ^  2  2  ^ 

i[p(l“P  )*csint  -(1-p  c  }cos0cost] 
2  2  ^ 

0  (1+  p  c  }cost  sin  0 


*  2  4 

-0[pccost  cos  0+(l-p  )*sin€] 


(22) 

The  functions  Gj(p),  ^2^^^  ^^(P)  used  in  the  time  derivative  of 

the  integrals  (15)  in  the  calculation  of  the  complete  field  for  the  fixed  buried 
couple.  Note  that  where  for  the  buried  force  the  initial  P  and  S  field  is  a  simple 
discontinuity,  i.  e. ,  a  step  function  of  displacement,  the  couple  will  be  associated 
with  Initial  delta  function  displacements. 
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Section  4;  Summary,  with  a  description  of  the  singular  parts  of  the  surface 
displacement  for  the  general  point  source. 

In  the  earlier  parts  of  this  paper,  we  have  discussed  the  details  of 
calculation  for  the  surface  displacements  produced  by  a  point  source.  The 
usual  problem  in  which  this  information  may  be  useful  is  the  one  of  recognising 
the  type  of  source  which  produces  a  given  displacement  field.  It  is  in  this 
context  that  knowledge  of  initial  displacements  is  of  use;  we  shall  therefore 
conclude  this  paper  by  giving  results  for  the  singular  parts  of  the  surface 
displacement  field  for  the  general  point  source. 

The  velocity  potential  for  the  general  point  source  of  order  {n+  1)  is  a 
homogeneous  function  of  degree  -(n  +  1)  in  space  and  time  variables;  it  is 
given,  for  y  >  0,  by  the  integrals 


i  F  V 


,3  n  Tr  IV^CAi  +  Bl+Ck] 
i-^)  B[c2t-ajX-a^z-py  ‘^"l  *^“3  » 


12  2  2  i 


with  a  =  (m‘‘-QrJ' -  p  =  {l-aj^  . 


The  functions 


and 


V,  = 


1  (C2-ajVj-a3V3-«V2) 


V,  = 


2  ■  (O2  -  Vj  -  «3V3  -  PV2) 
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Introduce  the  effect  of  steady  motion  for  t  >  0,  the  source  being  forced  to  move 
with  the  point  x  =  v^t,  y  =  V2t,  z  =  v^t  ,  The  prototype  of  these  potentials 
is  related  to  the  force  p(XJ[  +  YJ.  +  Z  k);  when  this  acts  at  the  origin  we  have 
a  first  order  source  for  which 

-  =  Xa^  +  Yor  +  Za^,  8n^  Ac^  =  Zp  -  Ya^,  =  Xa^  -  Za^  , 


and 

Sir^Cc^  =  Ya^  -  Xp  . 

Other  point  source  fields  may  be  obtained  by  differentiation  or  integration  of 
the  point  force  field.  The  general  form  given  arises  because  for  the  functional 
dependence  given,  we  always  have  the  identities 

Ot  ^Ct 

_d _ _L  Jl  _l  __i  _a. 

8x  ~  “  c^  *  dz  ~  c^  8t  ' 


for  a  compressional  field  which  moves  in  the  direction  of  positive  y 

_d. _ 2.  J_ 

8y  "  "c^  at  ' 

and  for  a  shear  field 

JL  JL 

8y  “  c^  8t 

For  the  initially  buried  source,  with  a  velocity  restricted  to  subsonic  values, 
the  surface  field  contains  three  circular  fronts  on  which  the  displacement  field 
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2  2  ^ 

or  Its  derivatives  are  singular,  these  being  associated  for  r  =  (Cjt  -h)^  with 

2  ^ 

the  arrival  of  the  primary  P  waves,  for  r  =  Cjt  -  h{l-m  )^/m  with  the  arrival 

of  the  refracted  P  waves  and  hence  of  the  head  waves  of  shear,  and  for 
2  2  - 

r  =(c^t  -  h  with  the  arrival  of  the  primary  S  waves. 

For  the  first  order  source  the  first  arrival  is  a  step-function  P  wave,  given 
from  the  integral  11a  with  an  amplitude 

4ir^p(m^-p^)*  F(«j,  a^)  Vj(aj,  a^)  [2p(l -p^)*X  +  (l“2p^)l] 

R(p) 

where  for  z>  0,  6  =  arctan  z/x,  this  is  evaluated  with  =  p  cos  6, 

2  2  ^ 

=  P  sin  0,  p  =  m(l-h  /c^t  )*  and  with  positive  values  for  the  radicals 
(1-p^)^  and  (m^-p^)*  . 

For  the  shear  field  there  are  the  distinct  situations  linked  with  the  ranges 
2  -  2 

h<  <  h/(l-m  )*  or  C2t>h/(l-m  )*  .  For  the  former  range,  the  refracted 
P  field  has  not  separated  itself  from  the  incident  S  field,  there  is  thus  no 
head  wave  of  shear,  and  the  arrival  of  shear  is  marked  by  a  step  function  of 
amplitude 

4it^(1-P^)*  “3)[Ccos  0-Asin0][p(l-2p^)X-p^(m^-p^)*J.] 

<^2 

4ir^  “3)  B  - 


9 


■64- 


#386 


where  V^,  A,  B  and  C  are  evaluated  at  the  points  =  pcos  0,  =  psln0 

2  2  2*^ 

with  p  =  (1-h  /c_t  .  When  head  waves  are  present,  there  is  no  discontinuity 

b  * 

in  the  displacement  field  associated  with  their  arrival,  this  for  the  first  order 
source.  The  same  step  function  is  present  at  the  arrival  of  the  primary  S  waves, 
but  it  is  clearly  masked  by  the  conjugate  logarithmic  singularity  associated 
with  the  head  waves.  This  contribution  has  the  form 

16tr  Aj[V2(Ccos  0  -  a  sin  0)]  {i(l"2p^)^  +  P(l“2p^)0“P^)*l}(P^“‘"^)* 
c^  [(l-2pV+  16p‘^(l-p^)(p^-m^)J 

2  2  2  - 

this  just  outside  the  circle  r  =  (c^t  -h  )^,  with  V,  A  and  C  again  calculated 

2  2  2 

at  the  points  =  pcos0,  =psin0,  when  ps^^l-h  /c^t  )*>m  and  with  the 
logarithmically  large  quantity  Aj  defined  by  the  equation 


Aj  =  arccosh 


0^1(1 -m^)^  -  h 
mCr^  +  h^-c^t^)^ 


These  are  the  singular  arrivals  for  the  first  order  source.  Each  differentiation 
increases  the  order  of  the  source  and  makes  the  singular  contributions  more 
singular.  For  the  source  of  order  (n-i- 1),  the  leading  singularity  linked  with 
the  arrival  of  the  head  wave  is  the  (n  -  l]^-th  derivative  of  a  step  function,  i.  e. , 
the  function  6^*^  with  an  amplitude  proportional  to  the  vector 


,  2. 

J[(l  -  2m  )  +  m  r 
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The  arrivals  associated  with  the  principal  wave  fronts  are  more  singular.  The 
primary  arrivals  are  dominated  by  the  nth  derivative  of  a  step  function  i.  e. , 
the  function  6^*^  while  the  logarithmic  head  wave  singularity  leads  to  the 
conjugate  function  of  degree  -n,  a  continuous  function  which  is  0(d  ”)  as 
the  distance  d  from  the  S  wave  approaches  zero  in  the  head  wave  region. 

The  relative  magnitudes  of  the  various  components  of  displacement  are  the  same 
as  for  the  first  order  source. 

If  the  source  velocity  is  supersonic,  we  expect  conical  singular  surfaces 

to  trail  behind  the  source  point  and  to  be  tangential  to  the  spherical  wave 

surfaces  linked  to  the  setting  up  of  the  source.  Only  if  the  source  point  is 

2  2 

outside  the  sphere  [r  +  (y-i-h)  =  c^t  and  above  the  tangent  surface 

2  2  Z  ^  2  2 

r(Cjt  -h  =  c^t  -yh  will  the  P  cone  cut  the  plane  y  =  0  to  show  a  new 

singular  wave  on  the  surface;  this  will  be  a  segment  of  an  ellipse  if  the  source 

point  is  below  the  lowest  point  y  =  -  h  -  Cjt  of  the  corresponding  spherical 

P  wave  (1.  e. ,  if  v^  <  -c^t)  otherwise  it  is  a  segment  of  a  hyperbola.  Likewise 

only  if  the  source  point  is  outside  the  spherical  S  wave  and  above  the  tangent 
2  2  2  ^  2  2  2 

surface  rCc^t  -h  )*  =  c^t  -h  will  the  S  cone  produce  its  effect  on  the  free 
surface,  and  again  if  ^2  <  “C^  there  is  a  singular  wave  of  elliptical  form  on 
the  surface,  and  this  becomes  hyberbollc  if  v_  >  -  c  . 

These  singularities  in  the  form  of  conic  sections  are  of  the  type  found  in  the 
description  of  steady  conical  fields.  For  the  P  wave  due  to  a  source  of  order  1  , 
we  find  specific  displacements  on  the  singular  front  in  the  form 
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-4ir  v^G{,ci^,a^)  ^2 


a  (c^Ct-l/a)  -(m  -c^/a  }(&  +n  )J* 


2  Z  ^  fc»  •- 

where  =  (v^  +  v^)*  ,  a  =  (v“+  v“)^  ,  and 


2.  2i 


2  2  *3) 


G(«j,  «3)  =  [2«(ajl  +«3]c)  +  (l-2«j  -2«3)l] 


with 


®1  ~ 


VlV^'l'|-VlV;»;-»V3'l'3 


av 


»  ®3  = 


av 


and  with 


(a^-cf)*; 

a*!*!  =  c^,  *|<3  =  Cj(ta-|)  * 


4* 


/  2  2,4 

-gj)  ^ 

2  -  Cj(ta-e) 


for 


and 


a^  »  VjX  +  v^h  +  V32  , 

2 

av^H  « -VjV2X+v^h-V3V2Z  , 

at  s  -V3X  +  VjZ 


The  singular  displacement  for  the  conical  S  wave  is  of  the  same  type.  We 
take  the  same  values  for  and  a^,  but  we  are  to  evaluate  the  expression 
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4n  v^H(a^,a3) 


a^[c^  (t  -  e/a)2  -(1  -  c^/a^)  (^2  ^  ^Z^ji 


when 


#  2  2,4 

(a  -03)^ 


/  2  2,4 

(»  -Oj)' 


*1=V"'  +3=c,Ut-i)  i-z  = ‘SuTrir  ” 


with 


R(®j,  “3)  = 


2(Ca^-Ag^)  {(g^i  ttt^k)[l-2Qr^-2a^-2ttP]  -!«> 
R(«l»  *3) 


(Ba3-CP)i+(Ap-Ba3)  k 

P 

2  2  2  4 

For  a  range  of  values  which  makes  the  radical  (m  -a^  -«3)*  Imaginary,  there 
is  a  head  wave  contribution:  this  fixes  the  head  wave  singularity  in  the  form 

[(1  -  Cj/a^X;^  +  D^)  -  Oj  (t  -  l/a)^]* 

The  arrival  of  the  head  wave,  linked  with  the  arrival  of  refracted  P  waves,  is 
not  marked  by  discontinuities  of  displacement. 

For  the  higher  ordor  source,  differentiation  with  respect  to  t  is  necessary; 
the  source  of  order  n+  1,  moving  steadily  with  supersonic  velocity,  has  conical 
wave  fronts  on  which  the  singularity  is  0(d"""^)  . 
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As  a  final  contribution  we  shall  discuss  the  source  which  appears  on  the 
surface.  In  addition  to  the  limiting  forms  of  expressions  which  are  found  for 
the  buried  source,  the  effect  of  the  singularity  due  to  the  Rayleigh  function  is 
to  give  extra  contributions  in  the  form 


u„  = 


8ir^{(l-2R^)Rl(FV^)l~  2R(R^-m^)^  Im(FV^)X } 
c^dCRj  (c^t^-rV 


and 


(2)  I6ir^  (R  [1~2R^+  2CR^-l)^(R^-m^)MRl(*|<) f  t  (R^~mV  Im(>ii)l 
“SR  „  nfu\  (cV-r^)2 


C2  D(R) 


within  the  circle  r  =  c^^t  .  Here, 

i}<  =  V2(C  cos  e  -  A  sin  0)  , 

is  evaluated  at  or^  =  R  cos  6,  =  R-  sin  6,  with  the  radicals  a,  p  taking 

positive  imaginary  values,  and  with  D(R}  defined  as  for  the  equations  (17)  . 

If  the  source  moves  on  the  surface  with  a  velocity  greater  than  c^^,  a 
wedge-shaped  Rayleigh  wave  is  set  up.  If  it  moves  with  a  supersonic  velocity 
there  are  conical  S  and  P  waves  set  up;  all  the  details  of  these  steadily 
moving  fields  may  be  picked  out  of  Part  2. 

The  surface  source  of  order  n-l- 1  sets  up  transient  P  and  S  waves  which 
are  0(d  ^  )  near  the  wave  front,  while  the  conical  P  and  S  waves  associated 
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—n 

with  steady  supersonic  motion  are  0(d  ^)  .  Thus  for  the  surface  source 

the  P  and  S  arrivals  are  of  the  same  singularity  as  the  head  wave  arrival 
for  the  buried  source. 

-n-- 

The  Rayleigh  wave  produced  by  the  setting  up  of  the  source  is  0(d 
near  the  wave  front,  this  being  comparable  with  the  arrival  of  steady  conical  P 
and  S  waves  for  a  buried  source.  The  strongest  singularity  is  that  of  the 
steady  wedge-shaped  Rayleigh  wave  for  the  moving  source,  this  being  0(d 

All  other  parts  of  the  surface  displacement  field  are  infinitely  differentiable. 
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Figure  la:  The  formal  integration  path  passing  below  singularities  at  p  -R, 
-1,  -m,  and  above  those  at  p  =  m,  1  and  R,  together  with  the  actual 
hyperbolic  path  of  integration,  deformed  to  pass  round  conjugate  complex 
zeros  of  the  wave  function.  Residues  at  these  conjugate  points  determine 
the  field  inside  the  characteristic  wave  envelope. 


Figure  lb:  When  the  locus  of  complex  zeros  passes  between  the  points  p  =  m, 
and  p  =  1,  the  actual  integration  path  must  be  composed  of  a  hyptfbola 
with  a  horizontal  loop  round  the  In’anch  point  p  =  m  .  The  complex  ztf  os 
of  the  wave  function  again  determine  the  field  inside  the  characteristic 
wave  envelope,  while  real  zeros  determine  the  field  in  a  head  wave  region. 
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Figure  3:  The  structure  of  the  field  due  to  a  primary  source  of  shear.  Two 
stages  are  shown.  Before  the  wave  front  of  the  primary  field  begins  to 
travel  at  the  critical  angle,  the  wave  fronts  both  of  reflected  S  and 
refracted  P  waves  travel  together  along  the  free  surface.  At  a  certain 
instant  the  refracted  P  wave  meets  the  free  surface  at  right  angles, 
and  it  may  only  continue  to  travel  at  the  velocity  of  dilatation  waves  by 
breaking  away  from  the  incident  and  reflected  S  waves.  Head  waves 
of  shear  now  appear  in  order  to  satisfy  boundary  conditions  at  the  surface. 
The  two  stages  shown  indicate  the  initial  development,  and  the  subsequent 
breakaway  of  the  refracted  field. 


